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ABSTRACT 

We present a detailed and unique mass budget for the high-surface-brightness galaxy UGC 463, 
showing it is dominated by dark matter (DM) at radii beyond one scale length (ha) and has a 
baryonic-to-DM mass ratio of approximately 1:3 within 4.2/i#. Assuming a constant scale height (h z , 
calculated via an empirical oblateness relation), we calculate dynamical disk mass surface densities 
from stellar kinematics, which provide vertical velocity dispersions after correcting for the shape of the 
stellar velocity ellipsoid (measured to have ag/aR = 1.04 ± 0.22 and <j z /<jr = 0.48 ± 0.09). We isolate 
the stellar mass surface density by accounting for all gas mass components and find an average K- 
band mass-to-light ratio of 0.22 ± 0.09(ran) ta^{sys) M G /£g ; iZibetti et all and |Bell_eLaD predict, 
respectively, 0.56 and 3.6 times our dynamical value based on stellar-population-synthesis modeling. 
The baryonic matter is submaximal by a factor of ~ 3 in mass and the baryonic-to-total circular-speed 
ratio is 0.6Y^_QQg(rari) lio'if (sj/s) a * 2.2ft.jj; however, the disk is globally stable with a multi-component 
stability that decreases asymptotically with radius to Q ~ 2. We directly calculate the circular speed 
of the DM halo by subtracting the baryonic contribution to the total circular speed; the result is 
equally well described by either a Navarro-Frenk- White halo or a pseudo-isothermal sphere. The 
volume density is dominated by DM at heights of \z\ > l.6h z for radii of R > /jr. As is shown 
in follow-up papers, UGC 463 is just one example among nearly all galaxies we have observed that 
contradict the hypothesis that high-surface-brightness spiral galaxies have maximal disks. 
Subject headings: dark matter — galaxies: fundamental parameters — galaxies: individual (UGC 463) 
— galaxies: kinematics and dynamics — galaxies: spiral — galaxies: structure 



1. INTRODUCTION 

A primary goal of modern extragalactic astronomy is to 
reduce the complex, stochastic process of galaxy forma- 
tion into a few fundamental physical parameters. Such 
a goal appears tractable given the tight scaling relations 
exhibited by galaxies over a large dynamic range in ob- 
served properties, which to fir st order may be tie d to a 
single physical characteristic (|Disnev et al.l 12008). For 
example, measures of galaxy size, luminosity, and a viri- 
alized dynamical quantity (such as the circular velocity in 
rotationally supported systems and velocity dispersion in 
pressure-dominated systems) demonstrate strong covari- 
ance. Correlations among galaxy prop erties are found 
in two-dimensional scatter plots (e.g.. iCourteau et al.l 
120071: INair et al.ll20Tol:[Saintonge fc Spekkensll201lD. lines 
throu gh multi-dimensional space (e.g.. iTollerud et al.l 
1201 11), and more complex , multi-dimensional manifolds 
(e.g. JZaritsky e t al. 2008). Empirical and theoretical un- 
ders tanding of these rela tions over cosmic time (as in, 
e.g., iDutton et aT1l2011aft are critical. 

Two long-sta nding scaling relation s are the Tully- 
Fisher relation (jTullv fc Fished 119771 hereafter the TF 
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relation) — the correlation between the rotation ve- 
locity of spiral galaxies and their total luminos i ty - 
and the Fundamental Pl ane (FP; iDressler et al.l 119871 : 
iDiorgovski fc Davislll987D — the plane relating the size, 
surface brightness, and velocity dispersion of ellipti- 
cal galaxies. These fundamental relations are strongly 
linked to mass: The baryonic TF (BTF) relation 
(|McGaugh et al.1 [2001 iMcGaugbl |2005| ) . created by re- 
placing total luminosity with total baryonic mass, ex- 
hibits less scatter than the nominal TF relation over a 
wide range of luminosity and spiral type. The mass- 
based FP ([Bolton et al.l l2007j) . incorporating the total 
(baryonic+dark- matter [DM]) mass surface density in- 
stead of surface brightness, also exhibits lower scatter 
than its luminosity-based counterpart. It is interest- 
ing that the residuals are reduced for both the BTF 
and mass-based FP relation despite the exclusion of 
DM from the former. The tightness of the BTF im- 
plies that either DM is a rather negligible mass com- 
ponent or there exists a strict proportionality, in both 
relative amplitude and distribution, between DM and 
baryonic mass in spiral galaxies. The former is incom- 
patible with our current understanding of gravity and 
the current pa radigm of hierarchical d i sk-galaxy forma- 
tion (see, e.g.. iFall fej Cfstathio ul 119801; [Dalcanton et all 
119971 : IMo et alJll998t lAgertz et alJl201lh . and the latter 
is tantamount to the discomforting disk-halo conspiracj0 

6 The observed fine tuning of the relative fraction and distri- 
bution of baryonic and DM mass required to produce a total ro- 
tation curve that is dominated by baryonic matter at small radii 
with a smoot h transition to a roughly constant rotation speed at 
large radii (cf. Cascrtano & van Gorkom 1991; Amorisco & Bcrtin 
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(|van Albada fc Sancisilll986l hereafter lvAS86) . One can 
begin to address this contentious issue by placing direct 
constraints on the detail mass composition of galaxies. 

Although there are multiple methods of measuring the 
total mass enclosed within a given radius (e.g., dynam- 
ics, lensing), a robust decomposition of total mass into 
fractional contributions from DM, stars, and the inter- 
stellar medium (ISM) is non-trivial. Measurements of 
the atomic-gas mass can be made directly using 21cm 
H I emission, and molecular-gas mass can be estimated 
using CO emission. However, stellar mass estimates 
depend on the calibration of stellar mass-to-light ra- 
tios, Y», via resolved stellar populations in the most 
nearby (dwarf) galaxies or stellar-population-synthesis 
(SPS) modeling of integrate d light. The latter re- 
mains substantially un certain (jConrov et al.l [2009. 2010; 
IConrov fc Gunnl[20Toh r 7 l 

Rotation-curve mass decompositions provide upper 
limits on T* when one adopts the "maximum-disk" 
hypothesis, the assumption that the rotation velocity 
at the center is domi nated by the l u minou s matter 
(|vAS86t) . For example, iBell fc de Jo"nel (|200ll) used the 
"maximum-disk" rotation c urve decompositions made by 
iVerheiienl (|200U hereafter I VP If ) to place limits on the 
allowed Y*. However, rotation curves cannot provide 
unique measurements of T„ as we hav e recentl y illus- 
trated (jBershadv et al.ll2010al hereafter iPaper J ): infer- 
ence of T* based on rotation-curve mass decomposi- 
tions are unconstrained due to the disk-halo degeneracy 
(|van Albada et al.lll985f ). 

Given the uncertainty in SPS model zero-points and 
the disk-halo degeneracy, a direct measurement of T* 
is nee de d. Following the work of Bahca ll fc Casertand 
(11984D. Ivan der Kruit fc Freeman! (11984 I19 86T). and 
iBottemal (|1993f). the DiskMass Survey (DMS; iPapeTI ) 
aims to tackle this problem via dynamical measurements 
of the mass surface density, Edyn, of >40 low-inclination, 
late-type galaxies. Our measurements uniquely describe 
the baryonic mass distributions and DM-halo density 
profiles, puM, of each galaxy within ^3 disk scale lengths 
Qir), thereby breaking the disk-halo degeneracy and al- 
lowing for detailed calculations of disk-galaxy mass bud- 
gets. In this paper, we focus on providing a detailed, 
initial example of these calculations using UGC 463, lo- 
cated at equatorial (J2000.0) coordinates (RA,DEC) = 
(00 h 43 m 32!5,+14 d 20 m 34 s ). We continue our series by 
summarizing the baryonic mass fractions in 29 additional 
galaxies in Bershady et al. (2011; hereafter Paper V), 
submitted. 

Here we summarize some salient properties of UGC 
463: It is a well isolated galaxy with a moder ately-hig h 
extrapolated central disk surface brightness ([Paper If ). 
which is a fac tor of ~ 2 above the mea n derive d by 
iFreemanl (|1970f ). It is of late type (SABc; IPaper J ) and 
demonstrates an interesting three-arm spiral structure. 
The SDSS g-band surface photometry demonstrates a 

2010). 



7 Here, the remarkable success of McGaugh (2005) in reducing 
the residuals in his BTF relation by usi ng Y* as derived from the 
mass-discrepancy— acceleration relation (McGaugh 2004) is note- 
worthy; however, it is possible t hat this is more reflective of the 
ability of MOND (Milgrorn 1983) to fit rotation curves and/or the 
disk-halo conspiracy than it is of the absolute calibration of these 
T* measurements. 



clear Type II s urface-brightness profile, as defined by 
IFreemanl (|1970f) . with a profile "break" at a radius of 
~ 15", well within the field-of-view (FOV) of our kine- 
matic data; the break becomes less pronounced toward 
longer wavelengths. The disk is also bright in the mid- 
and far-infrared Spitzer bands, suggestive of significant 
star-formation activity and molecular gas mass. In gen- 
eral, UGC 463 is unexceptional in its optical and near- 
infrared (NIR) color, size, and luminosity; however, it is 
slightly redder and more luminous (in Mk) than typical 
of galax ies in the DMS Phase-B sample (as defined in 
IPaper II ). 

Our study of UGC 463 is a detailed example in the 
use of our full suite of data to produce quantities of 
fundame ntal relevance to the science goals of the DMS 
(IPaper If), following mu ch of the f ormalism developed in 
iBershady et al.l (|2010bL hereafter IPaper III ) . Given the 
large number of observational ingredients, we have relo- 
cated some detailed information to future papers, which 
we refer to throughout our discussion. An outline of 
our paper is as follows: Section [2] presents all the data 
products. We derive the on-sky geometric projection of 
the disk using our two-dimensional kinematic data in 
Section 02 including an extensive discussion of the in- 
clination, i. Based on this projection geometry, we pro- 
duce azimuthally averaged kinematic profiles and beam- 
smearing corrections and discuss the axial symmetry of 
the galaxy in Section [4] In Section [5j we derive salient 
properties of the disk including the shape of the disk stel- 
lar velocity ellipsoid (SVE), the disk stability, mass sur- 
face densities of all baryonic components, and dynamical 
mass-to- light ratios. In Section [SJ we produce a detailed 
mass budget for UGC 463 out to 15 kpc (~ 4.2 scale 
lengths); this analysis relies on a traditional rotation- 
curve mass decomposition but uses our direct measure- 
ments of Y». Having established the mass distribution of 
all the baryonic components, Section [B] also presents the 
DM-halo density and enclosed-mass distribution. There- 
fore, Sections [2] -[4] are largely concerned with data han- 
dling, whereas Sections [5] - [6] produce the scientifically 
motivated calculations that result from these data. We 
summarize our study in Section [7j 

We note here the nomenclature e(x) signifies the mea- 
surement error in a;, a 7 is the azimuthal average of x, 
and (x) is the combined radial and azimuthal average of 
x. When quoting two sets of errors in any quantity (as 
done in the Abstract), the first and second set provide 
the random and systematic uncertainties, respectively. 

2. OBSERVATIONAL DATA 

The DMS has collected an extensi ve suite o f data to 
reach our science goals, as described in IPaper 1 We draw 
upon a large fraction of those observations specific to 
UGC 463 for use in this paper. Table CO provides a list 
of all the data products used herein, their observational 
source, the year of the relevant observations, a reference 
to the section containing a description of each dataset, 
and (when available) a reference to papers containing 
more detailed information. On-sky maps of much of the 
rele vant data products are provided in Figure CQ see Sec- 
tion 13.31 and Appendix [X] for a full description of how 
these maps were generated. 

2.1. Distance 
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Figure 1. Two-dimensional data used in our study of UGC 463. From top-to-bottom and left-to-right, we provide the surface brightness 
in AB mag arcsec -2 for SDSS g-band (n g ), SparsePak Ho-region (MHc*)' ant ^ SparsePak+PPak Mgl-region (^Mgi)i LOS velocity in kms — 1 
for Hi (Vjjj), Ho (Vrc), and stars (V*); mass surface density in Mq pc — 2 for Hi (Sjjj) and H2 (Sh 2 ); and stellar velocity dispersion 
((T*) in kms — 1 . The galaxy center, as provided by NED, is marked by either a black or white plus sign. All images have the same spatial 
scale and have a sky-right orientation. Contour levels are: fiua — 22, 21, 20 mag arcsec -2 ; /iMgi — 22.7, 21.7, 20.7 mag arcsec -2 ; Vjjj, 



4360, 4390, 4420, 4450, 4480, 4520, 4540 kms~ 



-Hi 



kms -1 . Gray contours are used to ease visibility with respect to the background 
(black) arc ovcrplotted on the fi g image for comparison. 



-1 



9 Mq pc -2 ; S H2 — 15, 30, 60 Mq pc -2 ; and a, — 30, 45, 60 
rayscale image. The contours of fm a (gray) and /%[gi 



The distance to UGC 463 is used to calculate: (1) 
the total absolute if-band magnitude, Mk, providing 
an inclination measurement via inversion of the TF re- 
lation (Section 13.1. 2[k and (2) the disk scale height 
based on a measure d scale length in kpc (see equa- 
tion 1 from IPaper ij ) . In Section 13.1.11 we find V sys = 
4460 ± 1 kms -1 , consistent with V sys = 4452 ± 9 kms" 1 
(jHuchra et al.l 119991 ) provid ed by NEDg| App lying the 
104 kms -1 flow correction (jMould et al J 12000(1 . we cal- 



8 The NASA/IPAC Extragalactic Database, operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Administration. 



culate a flow-corrected velocity of Vh w = 4356 ± 52 
kms -1 , where we have taken half the flow correction as 
its error ([Paper IIP . Using Ho — 73 ± 5 km s -1 Mpc" 1 
for Hubble's constant (pr ovided by NED, cf. iRiess et al.l 
l2009t lLarson et al.l 1201 If) , we calculate a flow-corrected 
distance of D = V aa w/H = 59.67 ± 0.01 ± 4.15 Mpc; the 
systematic error is dominated by the uncertainty in Hq. 

2.2. Optical and Near-Infrared Emission 

We use archival g- 1 r-, and i-band d ata obtained from 
the Sloan Digital Sky Survey (SDSS; lYork et alJl2000h 
and J-, H-, and if-band data ob tained from the Two- 
Micron All-Sky Survey (2MASS; ISkrutskie et all 120061) 
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Table 1 

Observational Data 



Product 



Source Year 



Section Ref a 



Ha kinematics 
[O III] kinematics 

stellar kinematics 

g, r, i photometry 
K, JHK photometry 
Hi aperture-synthesis 
24- (im photometry 



SparsePak 
SparsePak 
PPak 
SparsePak 
PPak 
SDSS 
2MASS 
VLA 
Spitzer 



2002 
2006 
2004 
2006 
2004 
2003 
2000 
2005 
2007 



2.3.1 


2.3.3 


2.3.2 


2.3.3 
233 


2.4.1 


2.4.3 


2.4.2 


2.4.3 




[2~2 

W2 
|2~5 


12.6.1 



a References: 1. S wa t ers et al., in prep; 2. Ande rsen et 
al., in prep; 3. IWestfalll lj2009l ); 4. iMartinssonl l|2011f ) 

to produce surface-brightness profiles and large-aperture 
total magnitudes. Photometric measurements are in AB 
magnitudes for SDSS data and Vega-based magnitudes 
for 2MASS data. SDSS and 2MASS images are, respec- 
tively, 10.'2 x 13^8 and 8f3 x 17.'1 with UGC 463 well 
separated from the frame edges. 

2.2.1. Surface Photometry 

Given the basic image reduction and photometric cali- 
bration provided by SDSS and 2MASS, our surface pho- 
tometry is primarily concerned with sky-background sub- 
traction and masking sources other than UGC 463. 

Source catalogs have been created for each band us- 
ing Source Extractor^] Each catalog has been visually 
inspected and pruned of erroneous source identifications, 
such as along meteor streaks or diffraction spikes; these 
features are masked from our final results by including 
pseudo-sources in our catalog. We have created a master 
catalog for the region surrounding UGC 463 by merging 
the catalogs from all bands, identifying sources detected 
in multiple bands. 

Using the IRA^\ task imsurf it, we determine the 
sky background of each image by fitting a Legendre poly- 
nomial surface (with cross terms) to each image where 
all sources and artifacts are replaced, initially, by a (±3) 
sigma-clipped mean of the image. After the lowest-order 
surface fit (2-x2-order), masked regions are replaced by 
the fitted surface values as the fit order is increased. By 
inspection, we find there is little improvement in the sky 
flatness when using surfaces of more than 9-x9-order 
(terms up to a: 8 ), and higher order fits begin to introduce 
artificial structure. The backgrounds of SDSS images are 
generally well-behaved, whereas the 2MASS H- and K- 
band data exhibit significant background structure. 

Surface photometry has been performed on each image 
after subtracting the sky background and masking all ar- 
tifacts and sources, except for UGC 463. Source mask- 
ing is forced to be identical in every band. From these 
masked images, we perform elliptical aperture photome- 
try over a range of radii, each with an aspect ratio and 
orientation coinciding with the derived geometry of the 
disk discussed in Section |3~T1 Given the shallow depth of 
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10 IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatory, which is operated 
by the Association of Universities for Research in Astronomy, Inc., 
under cooperative agreement with the National Science Founda- 
tion. 
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Figure 2. Optical and NIR photometry of UGC 463, corrected 
for Galactic extinction; SDSS measurements are in AB magnitudes 
and 2MASS measurements are in Vega magnitudes. Top — From 
top-to-bottom, the 2MASS K-, combined 2MASS JHK-, SDSS i-, 
SDSS r-, and SDSS g-band surface-brightness profiles are plotted as 
solid lines; the bands are differentiated by grayscale as given by the 
key. Exponential-disk fits with a fixed scale length of tiR = 12''3 
are shown as dashed lines (compare to the results in Table [2] when 
fiR is allowed to be free for each band); the vertical dotted lines 
mark s the innermost radius included in the fit. The A"-band data 
from dc Jons & van dcr Kruit (1994) are plotted as open circles. 
Bottom — SDSS g — r, r — i, and g — i colors; the thick gray 
line is the g — i color after correcting for internal extinction. The 
estimat e of t he internal i\-band extinction, A l j^ 27 , calculated in 
Section l2.2.4l is also provided as a dashed line. 

the 2MASS data, we have also produced a JHK surface- 
brightness profile using the unweighted sum of the J, 
H, and K surface-brightness profiles. This extends the 
NIR surface-brightness profile to larger radii. The 11 JHK 
bandpass" has an effective band-width of SX/X — 0.62 
and a Vega zero-point of 1062 Jy. 

Figure [2] provides all the surface photometry used in 
our present study of UGC 463. We apply Galactic extinc- 
tion correctionfOof A g = 0.357, A r = 6.242, A, = 0.173, 
and Ak = 0.034; no extinction correction is applied to 
(J- jhk (see Section |2"X21. We also plot the UGC 463 K - 
band photometry from Ide Jong fe van der Kruitl (fl99l . 
which is in very good agreement with our 2MASS pho- 
tometry at R < 40". Table [2] provides the result of fitting 
an exponential disk to all bands, including the JHK data, 
demonstrating marginal change in the best-fitting scale 
length and no evident trend with wavelength. The Ta- 
ble provides the radii over which the exponential disk is 
fit; the minimum radius is always 15" in order to avoid 
non-exponential features seen near the galaxy center. 

2.2.2. K-Band Surface Brightness Profile 

We apply two corrections to /ijhk and use it as the pri- 
mary NIR surface-brightness measurement for our study 
of UGC 463: (1) a color correction to produce a more 
accurate if-band surface-brightness profile at large ra- 
dius and (2) an instrumental-smoothing correction. The 

11 http: //irsa. ipac . caltech.edu/applications/DUST/ 
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Table 2 

Exponential Disk Scale Lengths and Scale Heights 



Band 


^max 

(arcsec) 


MO 

(mag arcsec -2 ) 


h R 
(arcsec) 


h R 
(kpc) 


WW 


(kpc) 


9 


70 


20.16 ±0.07 


12.0 ±0.5 


3.5 ±0.1 ±0.2 


8.1 ±0.3± 2 {° 


0.43 ± 0.03 


r 


70 


19.58 ±0.06 


12.2 ±0.4 


3.5 ±0.1 ±0.2 


8.1±0.2± 2 ;° 


0.44 ± 0.02 


% 


70 


19.29 ±0.05 


12.6 ±0.3 


3.6 ±0.1 ±0.3 


8.2±0.2± 2 ;2, 


0.44 ± 0.02 


K 


55 


16.66 ±0.09 


11.9 ±0.3 


3.4 ±0.1 ±0.2 


8.0±0.2± 2 ;° 


0.43 ± 0.02 


JHK 


80 


17.18 ±0.05 


12.3 ±0.3 


3.6 ±0.1 ±0.2 


8.1 ±0.2±l° 


0.44 ± 0.02 



+0.11 
-0.08 
+0.11 
-0.09 
+0.11 
-0.09 
+0.10 
-0.08 
+0.11 
-0.09 



Notes. Columns are: (1) photometric band; (2) maximum fitted radius, the minimum 
radius is always _R m i n = 15"; (3) best-fitting central surface brightness; (4) be st-fitting scale 
length; (5) scale length in kpc; (6) oblateness calculated using equation 1 from lPaper III ; and 
(7) scale height. 



color correction compares fix and [ijuk from Figure [21 
such that it includes the Galactic extinction correction 
to [J* jhk ■ We find hk -If JHK = -0.47 ±0.03 at R < 40" 
with a maximum deviation from the mean of 0.05 mag- 
nitudes, which is at most ~ e(Mif)- The NIR color gradi- 
ents are small over this radius, with A( J — H) < 0.2 and 
A(J — K) < 0.1, and roughly consistent with no color 
gradient to within the photometric error. 

Our instrumental-smoothing correction effectively per- 
forms a one-dimensional bulge-disk decomposition of 
\J-jhk- We assum e the c e ntral light concentration intrin- 
sically follows a iSersid (|1963[ ) profile. After first sub- 
tracting an exponential surface-brightness profile fitted 
to hjhk data between 10" < R < 33", we model the 
central light concentration within R < 8" by a Sersic 
profile convolved with a Gaussian kernel, assuming the 
latter is a good approximation for all instrumental ef- 
fects. We find a best-fitting Sersic index of n = 1.5 
and an effective (half-light) radius of R e = 1"7. Ran- 
dom errors in this modeling are approximated by the 
root-mean-square (RMS) difference between the mea- 
sured and modeled central light profile. We assume that 
the "intrinsic disk surface brightness" is the remainder 
of the profile after subtracting the model of the central- 
light concentration, the "intrinsic central light concen- 
tration" is the Sersic component of the model, and the 
"intrinsic \ijhk profile" is the sum of these two compo- 
nents; our instrumental-smoothing correction then con- 
sists of the difference between the measured and "intrin- 
sic" fijHK profile. We adopt a conservative 50% system- 
atic error in this correction due to the inherent uncer- 
tainties in the true parameterization of the central light 
concentration and instrumental-smoothing kernel. Al- 
though relevant to our assessment of any central mass 
concentration in our calculations of Sdyn (Section I5.3[) . 
T^ 1 * (Section 15. 6p . and the mass budget (Section [5]), 
our instrumental-smoothing correction is immaterial to 
the fundamental conclusions of our paper. 

Hereafter, we refer to the corrected JHK measurement 
as \j! k \ ha always refers to the measurement based on 
this profile unless otherwise stated. The fully corrected 
\x' K profile is shown in, e.g., Figure 01 

2.2.3. Total Apparent K-Band Magnitude 

In Section 13.1.21 we use Mk to calculate the inverse- 
TF inclina tion o f UGC 463 based on the TF relations 
derived by IV011 Therefore, we calculate the total ap- 



parent magnitude, tuk, using the 2 MAS S if-band data 
according to the procedure used by I VP 11 : We measure a 
roughly isophotal magnitude at = 21.5 mag arcsec -2 
(the rough surface-brightness limit occurring at R < 55") 
and extrapolate to infinity based on the fitted exponen- 
tial disk (h R = 12"3 ± 0"3). We find m K (R < 55") = 
9.45 ± 0.01 and a correction of -0.09 mag for the extrap- 
olation of the disk. Accounting for Galactic-extinction 
(A K = 0.034), we find m K = 9.32 ± 0.02. 

2.2.4. Internal Dust Extinction 

Left uncorrected, fi' K may overestimate our dynamical 
mass-to-light ratios (Section l5.6[) due to internal dust ex- 
tinction in if-ba nd, A l K . Using the dust-slab model from 
iTullv fc Fouauel (p85TPl and i = 27° ± 2° (Section EU, 
we calculate the function A l ^ Tl (R) as shown in Figure 
H which has a maximum of A 1 ^ 27 (R = 0) = 0.02 ± 0.01 
at the galaxy cente r. Assuming Av/E{B — V) = 3.1 
(jCardelli et al.lll989D . we calculate A K /E(g-i) = 0.18 ± 
0.02, which we use to apply an internal reddening cor- 
rection to our g — i color when comparing our dynam- 
ical T* measurements to those predicted by SPS mod- 
eling. Figure [2] plots the uncorrected g — r, r — i, and 
g — i colors, as well as our internal-reddening-corrected 
(<7 — i)o color. Compared to more realistic radiative- 
transfer modeling including a clumpy ISM and spiral 
structure, these simple dust-slab model predictions tend 
to overestimate both the level of dust extinction and 
reddening toward high inclination; the predictions are 
more reasonable at the low inclination appropriate for 
UGC 463 (Schechtman-Rook et al., in prep). Given the 
marginal extinction in K-hand and reddening of g — i, 
the simpler model is sufficient for our purposes. 

2.2.5. Dust Emission 

We note here that the if-band also contains emission 
from hot dust and polycyclic aromatic hydrocarbons; 
however, based on a preliminary modeling of the spectral 
energy distribution of our full suite of NIR and Spitzer 
imaging, we expect this to be no more than a 3% con- 
tribution, which is immaterial to the conclusions of this 
paper. 

2.3. Ionized-Gas Kinematics 

12 See also discussion in lVOll and lPaper III . 
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We primarily use io nized-gas kinem a tics o b tained 
by the Spa rsePajEj (IBershadv eFall [200l 120051) 
and PPalgl (IVerheiien et al.l l200i IRoth et all [2001 
iKelz et al.f l2006D integral-field units (IFUs), augmented 
by our H I observations, to produce the total rotation 
curve of UGC 463. These ionized-gas data also provide 
measurements of the gas velocity dispersion, which we 
use to correct the gas rotation speed to the circular speed 
(Section 14.3. Below, we briefly describe the IFU data 
available for UGC 463 and our extraction of the ionized- 
gas kinematics. 

2.3.1. Ha, Nil, & Sll Spectroscopy 

SparsePak integral-field spectroscopy (IFS) of UGC 
463 was obtained on the nights of UT 02 January 2002 
and UT 20 October 2002, following the s etup pro vided 
for the Ha region as listed in Table 1 of IPaper 1 We 
obtained four pointings during the January run and an 
additional three pointings during the October run. The 
pointings nominally followed the 3-pointing dither pat- 
tern design ed to fully sample t he 72" x 71" FOV of 
SparsePak ([Bershadv et al.l I2004T )P1 the fourth pointing 
during the January 2002 run was a repeat of the cen- 
ter pointing. We obtained 2xl5-minute exposures for 
each pointing. Each exposure pair is combined, before 
extraction of the spectra, while simultaneously remov- 
ing cosmic rays. Spatially overlapping fibers among the 
seven pointings are not combined but treated individu- 
ally throughout our analysis. Further details of the re- 
duction of these data (basic image reduction, spectral 
extraction, wavelength calibration, and sky and contin- 
uum subtraction) are provided by S wate rs et al., in prep, 
largely following methods described in lAndersen et al.l 
( 2006) with continuum-s ubtraction techniques described 
in IBershadv et al.l ()2005[ ). The RMS difference between 
the catalogued and measured line centroids for the ThAr 
lines used in our wavelength calibration, i.e. the "wave- 
length calibration error," is typically < 0.1 kms -1 . 

We also measure the instrumental dispersion, (Ti ns t, 
as a function of wavelength for all spectra, using the 
ThAr em ission lines from our calibration lamp spectra 
(|Paper III ). The intrinsic widths of the ThAr features 
are negligible such that the second moment of these 
lines is equivalent to Cinst to good approximation. Af- 
ter identifying a set of appropriate (unblended) lines 
from the calibration spectrum, we fit single Gaussian 
functions to each line using the same code described by 
lAndersen et al. (2008) to fit the emission- line features 
in our galaxy spectra (see also Swaters et al., in prep). 
For each fiber, we fit a quadratic Legendre polynomial to 
Cinst(A), which is used to interpolate <7i, ls t at any wave- 
length. The average instrumental broadening across the 
full spectral range for all fibers is 13 kms -1 . 

2.3.2. Om Spectroscopy 

13 Mounted on the 3.5-meter WIYN telescope, a joint facility 
of the University of Wisconsin-Madison, Indiana University, Yale 
University, and the National Optical Astronomy Observatories. 

14 Mounted with PMAS on the 3.5-meter telescope at the Calar 
Alto Observatory, operated jointly by the Max-Planck-Institut fr 
Astronomic (MPIA) in Heidelberg, Germany, and the Instituto de 
Astrofsica de Andaluca (CSIC) in Granada, Spain. 

15 http: //www. astro .wise . edu/~mab/research/sparsepak/ 



Our optical continuum spectra in the Mg I region 
taken with both SparsePak and PPak — described in 
Sections 12.4.11 and 12.4.21 respectively — have sufficient 
spectral range to include the [O iii]A5007 emission fea- 
ture. Therefore, we also use these lines as tracers of the 
ionized-gas kinematics. No adjustment of the continuum- 
data reduction recipe was needed to accommodate the 
proper handling of the emission features. Instrumental 
dispersions are calculated as described in Sections 12.4.11 
and l2.4.2"l for the SparsePak and PPak data, respectively. 

2.3.3. Kinematic Measurements 

Ionized-gas kinematics are measured f o r all available 
emission lines. Follow ing lAndersen et al.l (|2008l see also 
lAndersen et~aT]|2006D . both single and double Gaussian 
line profiles are fitted in a 20A window centered around 
each line. All Gaussian fits have been visually inspected 
to ensure each emission line was fitted properly. Veloci- 
ties (cz) of each atomic species are calculated separately 
using the wavelength of the Gaussian centroid. Of all fit- 
ted line profiles, 27% are bett er fit by a double Gaussian 
profile (Ander sen et al.l [20081 Andersen et al., in prep); 
in these cases, a single component is used to measure the 
line-of-sight (LOS) velocity. Ionized-gas velocity disper- 
sions, Cgas, also use a single component and are corrected 
for the measured instrumental line width. 

For our Ha-region spectroscopy, we combine all avail- 
able velocity and velocity dispersion measurements 
(any combination of the [N n]A6548, Ha, [N n]A6583, 
[Sll]A6716, and [Sn]A6731 lines) into an error-weighted 
mean velocity for each fiber. Due to the large uncer- 
tainties in cTgas for lines other than Ha, we only include 
measurements with e(fr gas ) < 3 kms -1 in the combined 
value. In this spectral region, both velocities and a gas 
are dominated by the Ha line measurements due to the 
higher S/N of these lines; hereafter, we refer to these 
kinematics as "Ha" kinematics, despite their inclusion 
of other ionized atomic species. 

In addition to the consistency check among measure- 
ments made by SparsePak and PPak, the [O in] kine- 
matics provide a useful comparison with the Ha results 
(see Section [4}. However, the Ha line generally provides 
higher quality kinematics and velocity fields: the line 
S/N is higher on average and the filling factor of the 
kinematic measurements in the disk of UGC 463 is more 
uniform. We eventually combine all ionized-gas kine- 
matics into a single, axisymmetric set of measurements 
(Section I4.3|) ; however, it is useful to keep in mind this 
distinction between the merit of the Ha and [O in] kine- 
matic data. 

2.4. Stellar Kinematics 

Our optical continuum spectra are at the heart of this 
paper and, in fact, our entire survey. UGC 463 is among 
a set of 19 galaxies in (and roughly half of) our Phase-B 
sample that have both SparsePak and PPak continuum 
spectroscopy near the Mg I triplet. This intentional du- 
plication provides an internal consistency check of our 
stellar kinematics using two different instruments, and 
we find excellent agreement among the observations (Sec- 
tion Here, we provide information concerning our 
observations and our derivation of the LOS stellar kine- 
matics. 
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2.4.1. SparsePak Spectroscopy 

SparsePak IFS of UGC 463 was obt ained us ing the 
Mg i-region setup as listed in Table 1 of IPaper J — pri- 
marily targeting Fei and Mgi stellar- atmospheric absorp- 
tion lines. The o bservation and reduction of these data 
arc described bv IWestfalll (2009); we review the salient 
details here. 

UGC 463 was observed on consecutive nights during 
a single run from UT 23 — 25 September 2006. Four, 
six, and four 45-minute exposures were taken during the 
three nights of observation, respectively. No dithering of 
the pointing was applied between exposures; the original 
pointing was repeated to the best of our ability for each 
night. Exposures taken within a given night have been 
combined into a single image, using an algorithm that 
simultaneously rejects cosmic rays, and reduced (basic 
image reduction, spectral extraction, wavelength calibra- 
tion, and sky-subtraction) on a night-by-night basis. The 
basic reduction procedures are nearly the same as that 
used for the Ha data (Section 12.3. 1[) . The wavelength 
calibration errors are typically < 0.1 kms -1 . Error spec- 
tra have been calculated in a robust and parallel analysis. 
The repeatability of the pointing across the three nights 
of observ ation was good to less than one arcsecond, deter- 
mined by Westfall (2009) by forcing the kinematic centers 
of all the SparsePak Ha and [O in] data to follow from the 
same on-sky kinematic geometry. Thus, the extracted 
spectra from each night have been combined on a fiber- 
by-fiber basis and weighted by the spectral (S/N) 2 . The 
weighting is relevant due to changing conditions; moon 
illumination increased for each night (from 1 — 7%) and 
a number of exposures (~ 30%) suffered from variable 
transparency losses due to passing clouds, particularly 
during the second night of observation. The combined 
spectra, from all 10.5 hours of integration, are analyzed 
in Section 12.4.31 to measure stellar kinematics for UGC 
463. 

Using our calculated error spectra, Figure [3] plots the 
mean S/N of our SparsePak IFS as a function of radius 
against the g-band surface-brightness profile, demon- 
strating that we have measured stellar kinematics for 
fibers with fi g ~ 22.5 at S/N ~ 3. The scale transla- 
tion between mean S/N and assumes the data are 
detector-limited (S/N ex flux), which is only an approx- 
imation for our data. 

Finally, we estimate <Tj nst for both the template 
and galaxy spectra for use in measuring instrumental- 
broadening corrections (Section l2.4.4j) . We use the same 
approach as described for the Ha-region spectroscopy in 
Section 12.3.11 for both the galaxy and template obser- 
vations. The average instrumental broadening for these 
galaxy spectra is 11 kms -1 . 

2.4.2. PPak Spectroscopy 

PPak IFS of UGC 463 was obtai ned usin g the Mg I- 
region setup as listed in Table 1 of IPaper H Compared 
to the SparsePak optical continuum spectra, individual 
PPak spectra have approximately 1.5 times the spectral 
range, 0.7 times the spectral resolution, and 0.6 times 
the on-sky aperture per fiber; however, PPak contains 
4.4 times as many fibers as SparseP ak in its ma i n fibe r 
bundle. These data are described by Martinsson (2011), 
which includes all Mgi- region data taken by PPak. Here, 
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Figure 3. Mean S/N measurements for both our SparsePak (top) 
and PPak (bottom) Mgl-region IFS of UGC 463. Data are plotted 
as a function of rad ius and against the SDSS g-band profile (gray 
line; Section 12.2,11) , Bars provide the range in S/N of all pixels 
over the full spectrum; upper limit arrows are used in cases where 
S/N drops below zero due to the random errors in sky subtraction. 
Spectra that have yielded stellar kinematics are plotted by filled 
points and as open circles otherwise. 



wc briefly review the acquisition and reduction of the 
data specifically for UGC 463. 

We obtained two consecutive 1-hour exposures target- 
ing UGC 463 on each of three nights from UT 13 — 15 
November 2004, yielding six hours of total on-targct 
integration. Flexure corrections have been applied on 
an exposure-by-exposure basis, requiring spectral extrac- 
tion, wavelength calibration, and sky subtraction to be 
performed on each exposure individually. The RMS 
wavelength calibration error is typically of the same or- 
der as that found for the SparsePak observations (~ 0.1 
kms -1 ). Error spectra have been calculated in a robust 
and parallel analysis. For UGC 463, no pointing offsets 
were detected; the CCD of the guide camera of the PMAS 
spectrograph allows for accurate reacquisition of our tar- 
get galaxies on subsequent nights. All six spectra for a 
given fiber have been combined using weights depending 
on the instr umental reso l ution and the spectral S/N as 
described bv IMartinssonl f|201 ID - 

Figure [3] provides the mean S/N of the PPak spec- 
tra, alongside that of the SparsePak spectra. Despite 
the shorter integration time and smaller fiber aperture, 
the PPak data have slightly higher mean S/N, an effect 
of both the lower spectral resolution and the better effi- 
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ciency of PMAS over the WIYN Bench Spectrograph^! 

Instrumental broadening measurements are performed 
differently for PPak data than described above for 
SparsePak data. PPak observations provide simultane- 
ous calibration spectra in 15 fibers, evenly distributed 
along the pseudo-slit among the galaxy spectra. These 
spectra have prove n critical for apply ing the necessary 
flexure corrections ([Martinssonl 1201 if ) . Moreover, they 
provide simultaneous measurements of the instrumen- 
tal dispersion, obtained by fitting Gaus s ian fu nctions to 
the ThAr emission lines. iMartinssonl ([20111 ) has fit a 
quadratic Legendre polynomial surface to these measure- 
ments of CTi ns t for each object frame. We use this descrip- 
tion of (Ti ns t to calculate the instrumental-dispersion cor- 
rections for both the ionized-gas and stellar kinematics 
measured from the PPak spectra. The mean instrumen- 
tal dispersion across all fibers and all spectral channels 
is 17 kins" 1 for our PPak data of UGC 463. 

2.4.3. Raw Kinematic Measurements 

We use the Detector- C ensored Cross-Corre lation 
{DC S) software presented in iWestfall et~al~l (|201lL here- 
after iPaper HI ) to determine spatially resolved stel- 
lar kinematics in UGC 463 from both our SparsePak 
and PPak spectra. Based on a preliminary analysis 
(|Paper III ), we find a Kl III star provides a minimum 
template-mismatch error of ~ 5% in the observed veloc- 
ity dispersion (er bs) with no systematic trend in radius. 
In the future, we can improve upon this by using compos- 
ite templates; however, a 5% template-mismatch error in 
(T bs is satisfactory for the pr esent stu dy. 

As per our survey protocol (|Paper If ), we have observed 
template stars in the Mg I region using both SparsePak 
and PPak. For the analysis here, we specifically use 
the Kl III stars HD 167042 and HD 162555 for our 
SparsePak and PPak stellar kinematics, respectively; Ta- 
ble [3] presents salient information regarding the template 
spectra. Template star observations are performed un- 
der nominally the same spectrograph/telescope configu- 
ration as for our galaxy data. For both SparsePak and 
PPak, template stars are observed by drifting the star 
through the full FOV, yielding many spectra that are 
combined to provide high-S/N templates; the final S/N 
of each template is provided in Table |3l 

We fit any galaxy-template cross-correlation (CC) 
function that peaks within a few hundred kms -1 of the 
systemic velocity of UGC 463 (as recorded by NED) re- 
gardless of the S/N. For each fiber, we adopt a Gaus- 
sian broadening function, and we use a cubic Legendre 
polynomial to minimize continuum differences between 
the broadened template and the fitted galaxy spectrum. 
Additionally, we mask the [O m]A5007 and [N i] (A5198 
and A5200) nebular emission regions from both the tem- 
plate and galaxy spectra. For the PPak data, the red- 
shifted [O lil]A4960 line is also visible; however, DCS 
masks the CC to a rest-wavelength range common to 
both the galaxy and template spectrum, thereby auto- 
matically masking this line. Each CC fit has been visu- 
ally inspected to insure the proper peak was considered 
by the fitting algorithm and that any unexpected arti- 

16 Our data were taken before the completion of the upgrade to 
the WIYN Bench Spectrograph that improved its over all efficiency 
by a factor of > 2 ( Bersh adveFaT| [2008 : Kncz ek et aljEaoTOT) . 



facts — poorly removed sky lines and/or cosmic-ray de- 
tections — were masked. Based on this inspection, spec- 
tra have been refit as necessary. Our stellar kinematic 
analysi s follows t he expectations derived for random er- 
rors in IPaper III As assessed via \ 2 an d the velocity 
shift with respect to spatially neighboring fibers, we find 
reasonable fits to spectra with mean S/N approaching 
unity, albeit with large errors. Systematic errors should 
be negligible for velocity measurements at all S/N, and 
they should be < 20% in <7 b s at S/N > 2; systematic er- 
rors are alw ays smaller than the calculated random error 
(|Paper IIS ). 

2.4.4. Instrumental-broadening Corrections 

We correct our observed stellar kinematics for the sys- 
tem response function by considering the following two 
separable components: (1) The broadening of the intrin- 
sic absorption-line widths due to the spectrograph op- 
tics, accounted for using an "instrumental-broadening" 
correction, 5cr inst ; and (2) The smearing of the intrinsic 
surface-brightness, velocity, and velocity dispersion dis- 
tributions by the response of the atmosphere+telescope 
system, accounted for using a "beam-smearing" correc- 
tion, CTbcam- The final LOS dispersion is crj = a^ hs — 
Saf nst - crg cam ([Paper IS ). Unlike a bo am, Sa inst is inde- 
pendent of the on-sky geometry and intrinsic kinematic 
structure of the observed galaxy; therefore, we calculate 
^Uinst here. We calculate <7beam before combining our 
SparsePak and PPak kinematics in Section 2] using the 
projection geometry derived in Section O 

Each CC is used to compare template and galaxy 
absorption-line shapes such that <5cr inst is determined by 
the difference in <7i ns t measured for the template and 
gala xy spectru m; we calculate forest following Appendix 
A of IPaper "TP] using our measurements of CTi ns t for both 
the templ ate and g alaxy spectra. We adopt a 4% error 
in (5(7i ns t (jPaper IS ), which is marginal when compared 
to e(er bs). These corrections differ rather dramatically 
between SparsePak and PPak; however, in both cases, 
Saf nsi is typically small. Corrections to er b s — i-e., the 

ratio (al hs -Saf nst )^ /a obs — are < 4% and < 20% for, re- 
spectively, 90% and 99% of all measurements; a few mea- 
surements have rather large corrections due to dispersion 
measuremen ts of <7 bs < 10 kms -1 , which are likely erro- 
neously low ([Paper IIS ). We always find <5eri nst < e(cTi nst ). 

2.5. Atomic-Gas Content 

As part of our general survey strategy ([Paper J ), we 
have obtained 21cm aperture-synthesis imaging for the 
DMS Phase-B sample. These data measure neutral hy- 
drogen (H i) surface densities (Shi) an d extend the 
rotation-curve measurements of each galaxy; the ionized- 
gas kinematics can be limited by the FOV of our Ha 
spectroscopy and/or the extent of the Ha emission in 
the disk. For UGC 463, we obtained 2.3 hours of on- 
source integration using the Very Large Array (VLA); 
observations were taken in the C configuration yielding 
a synthesized beam of 14'.'7 x 12'.'9 and a velocity res- 
olution of 10.5 kms -1 . In the end, these data provide 
only a marginal radial extension of the ionized-gas rota- 
tion curve of UGC 463. The acq uisition and reduc tion 
of these data is fully described by iMartinsso n1 (j201~l . 

The two-dimensional Hi mass-surface-density map and 
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Table 3 

Stellar Templates 











S/N 
(pix-l) 


Physical Quantities 51 




HD 


Spectral 
Type 


Instrument 


UT Date 


Vhol '^eff log g 
(kms" 1 ) (K) (cms" 2 ) 


[Fe/Ej 
(dex) 


167042 
162555 


Kl III 
Kl III 


SparsePak 
PPak 


2001-06-09 
2007-01-15 


0.4xl0 3 
l.lxlO 3 


-18.01 ±0.17 4878 2.74 
-14.84 ±0.20 4660 2.72 


-0.11 
-0.21 



Measurements of V^el are from]Famacy ct al. (2005); remaining data are from lValdes et all 1120041 ). 
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Figure 4. Mass surface density measurements for Ht (via direct 
observation; open circles) and H2 (via an indirect calculation; gray- 
filled points), and all gas including helium and metals (crosses). 
The gray-filled points also provide the /24^m profile as delineated 
by the ordinate to the far left. Errors (random in black and sys- 
tematic in gray) shown for Ejj 2 and S gas are discussed in the text; 
errors in Sjjj are of order the size of the plotted symbol. As delin- 
eated by the right ordinate, we overplot [i' K (Section 12. 2.21 1 as the 
solid line with the error region plotted in dark gray. 



velocity field are presented in Figure [T] The azimuthally 
averaged measurements of Shi (R) f° r UGC 463 are pre- 
sented in Figure [4] we adopt c(£hi) = As i s 
typical of late-type spiral galaxies we find a decrease in 
the H 1 mass surface density toward the galaxy center; 
the peak surface density of Shi — 5.6 Mq pc -2 occurs 
at R = 1.9fiR. Although the spatial resolution of our 
H 1 column-density map is a factor of > 3 larger than 
our optical IFU data, we do not attempt to match the 
resolution of these two data sets; such a correction to 
the azimuthally averaged total mass surface density is 
negligible for the purposes of this paper. 

2.6. Molecular-Gas Content 

Our nominal survey scope does not include direct ob- 
servations of the molecular content of our galaxy sam- 
ple via, e.g., the 12 CO (J = 1 — > 0) emission line - 
henceforth all discussion of "CO emission" refers to this 
emission feature unless noted otherwise. Until we ob- 
tain such data, estimation of the molecular content in 
DMS galaxies relies on available literature data and/or 
inference from other suitable tracers for which data is 
available. 

For our characterization of the molecular content of 
UGC 463, we use our 24/im Spitzer imaging to produce 
a rough approximation of the CO surface-brightness dis- 
tribution. Numerous studies exist demonstrating a cor- 
relation between the infrared luminosity of star-forming 
galaxies (dominated by thermal dust emission) and their 
molecular-gas content as traced by CO emission. For 
example, the integrated infrared luminosity based on 



IRAS observation ^ 17 ! is well-correlated w ith the inte- 
grated CO fl ux (|Young fc Scoviflel 1199 If ) - Moreover, 
IBendo et al.l ((2007h note a similar dependence of the 
distribution of 24 um emission on m orphological type 
as was noted by lYoung et al.l (|1995f ) for the molecu- 
lar gas traced by CO emiss ion. At spatially resolved 
scales, iPaladino et al.l ()2006|) have studied the correla- 
tion between CO (IcoAV) and 24/im emission (/24/im) 
in a set of 6 nearby spiral galaxies to find /24/im oc 
(7coAU) - 9±01 . Correlations betwee n the CO and 8[im 
emission have als o been discussed (jRegan et al.l 120061 : 
IBendo et al.l I20l"0h ; however, we prefer to focus on the 
correlation between IqoAV and /24/im as the latter 
should be dominated by warm dust emission and be 
less dependent on the fraction of dust in the form of 
polyc yclic aromatic hydrocarbons (see, e.g-. IDraine k Li 
1200% . 

In Section 12.6.11 we describe the procedure used 
to measure /24/im (-R) in UGC 463. In Section 12.6.21 
we detail the calibration of our /24/tm-to-/co AU 
sur face-brightn ess relation usin g da ta made available 
bv ILerov et"ail (|2008l hereafter |L0l). Finally, in Sec- 
tion [5^31 we convert from /24/im to Ico^-V and then 
calculate the molecular mass surface density, £h 2 , using 
the traditional A-factor, Aco- This latter step domi- 
nates the systematic error in our estimation of the molec- 
ular content of UGC 463. 

2.6.1. 24fJ-Tn Spitzer Photometry 

The survey strategy for all our Spit zer observations are 
provided in Section 6.2.3 of iPaperll . In general, 24^m 
images collected for the DMS demonstrate significant 
background structure, due to both detector effects and 
intrinsic structure in the Galactic ISM, with fluctuations 
on angular scales close to that of our galaxies. To ac- 
count for these fluctuations, we mask out all statistically 
significant sources, including a substantial radial region 
surrounding UGC 463, and create a 47" x 47" boxcar- 
smoothed background image. Masked regions are itera- 
tively filled by the boxcar smoothing, effectively interpo- 
lating the sky background and its gross structure, across 
all detected sources. We simply subtract this smoothed 
image from our 24/itm image of UGC 463 and use the 
result to calculate the 24//m surface-brightness profile. 

Our background-subtraction procedure has been care- 
fully assessed to ensure that the low-surface-brightness 
extent of UGC 463 has not been systematically over- 
subtracted. Preliminary tests with UGC 463 and other 
galaxies in our survey demonstrate that our /24/im pro- 
files become strongly affected by the sky-subtraction er- 
rors at a source intensity below /24/im < 0-05 MJy sr _1 

17 http : //irsa. ipac . caltech.edu/IRASdocs/iras . html 



10 



Westfall et al. 



(Eh 2 < 0.21 M Q pc" 2 in Figure E]). UGC 463 is the 
third brightest 24/im emitter in our entire sample, mean- 
ing that this surface brightness limit falls outside the ra- 
dial region relevant to this paper. Our measured 24/im 
surface-brightness profile uses elliptical apertures follow- 
ing a geometry identical to that used for the optical and 
NIR photometry in Section 12.2.11 Figure 2] provides the 
24/im surface-brightness profile and the result of its con- 
version to Eh 2 , according to the discussion in the next 
two sections. The random errors in our /24^m mea- 
surements incorporate a constant 4% calibration error 
(jEngelbracht et aLll2007| ) and a sky-subtraction error es- 
timated by the change in /24^m introduced by a factor of 
two change in the smoothing-box size; the latter results 
in 1% and 10% sky-subtraction errors at R ~ 30" and 
~ 47", respectively. 

2.6.2. 24^m-to-CO Surface Brightness Calibration 

We use m easur ements of both CO and 24/im emission 
provided by IL08I (see their Table 7) to measure the cor- 
relation between J2 4 um and Ico&V- Twelve of the 23 
galaxies studied by IL08I include both CO and 24/im ob- 
servations; however, four of those galaxies (NGC 2841, 
NGC 3627, NGC 4736, and NGC 5194) are listed in 
NED as having either LINER or Seyfert activity, un- 
like UGC 463. Therefore, we calibrate /24/jm/JcoAU 
using only the remaining eight galaxies, hereafter the 
"h^m/Ico&V subsample." 

The quantities provided by IL08I are matched- 
resolution, azimuthally averaged radial profiles of £h 2 
(based on CO emission and a value for Xqq) and the 
contribution of embedded star formation (determined 
from the 24/im surface brightness) to the total star- 
formation-rate surface density. We revert these quan- 
tities to Ico^V (in K kms -1 ) a nd Io 4, im (in MJy sr _1 ) 
using equations A2 and Dl from IL081 All eight galaxies 
in the I^ian/ Ico&V subsa mple were observe d as part 
of the HERACLES Survey (ILerov et al.ll2009h . observ- 
ing only the 12 CO(J = 2 — > 1) emission line, where lL08l 
adopt a line ratio of 12 CO( J = 2 -> l)/ 12 CO( J = 1 -> 0) 
= 0.8. The CO surface brightness has been determined 
by integrating the emission profile over the full line width 
and converting the flux units per beam to Kelvin using 
the Rayleigh- Jeans limit. The 24/im fluxes are deter- 
mined from surface photometry of Spitzer imaging da ta 
obtained by the SINGS Survey (jKennicutt et al.l l2003). 

Figure [5] presents the data for the hA^m/Ico sub- 
sample regardless of the galaxy or radial regio n from 
which it has been measured. Table 7 from IL08I is used 
to calculate e(icoAU) directly, whereas we adopt a uni- 
form e(/24^m) = 0.15/24^im due to insufficient informa- 
tion; we expect this e(l24p m ) to be an upper limit. We 
fit a power-law relationship between and IcoAV 

to all available data, in corporating error s in both coor- 
dinates (Section 15.3 of lPress et al.ll2007t ) . finding a best 
fit of 



log 



Ico&V 
K km s" 1 



= 1.08 log 



'24^m 



MJy sr 



0.15 



(1) 

with a weighted standard deviation of ±0. 11 dex, in good 
agreem ent with the previous result from iPaladino et al.l 
(2006]). Thus, given that UGC 463 has physical param- 
eters that are comparable to the Iz^m/IcoAV subsam- 
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Figure 5. Correlation between azimu thall y averaged values of the 
24/im and CO surface brightness from L08 for the -/24u m /^coAV 
subsample (see text). The best-fit linear correlation (equation fTJ 
is given by the solid line; the dashed lines illustrate the weighted 
standard de viation about the fit. T he dotted line has the slope de- 
termined by Paladino et al. (2006) using a different galaxy sample 
with different CO data. The arrow at the top of the plot indicates 
the peak /24/jm measured for UGC 463. 

pie, our calibration is expected to estimate TcoAU for 
this galaxy to within ~ 30%. 

2.6.3. H2 Mass Surface Density 

Using equation [TJ we convert our sky-subtracted 24/im 
image of UGC 463 to a CO surface brightness map. Sub- 
sequently, we calculate the H2 mass surface density using 
the traditional conversion factor, Xco: following 



M Q pc" 



= 1.6 



IcoAV 



K km s 1 
Aco 



10 20 cm~ 2 (K km s" 



cos i, (2) 



where i is the galaxy inclination and Aco is the ratio of 
the H2 column density to the CO line strength. This use 
of Xqo to calculate Sjj 2 is a common procedure for cal- 
culating the molecular-gas content of external galaxies; 
however, it may suffer from substantial systematic error. 

A large number of studies have been devoted to mea- 
suring Xqo both in our own Galaxy and within the 
Local Group. Empirical and theoretical studies sug- 
gest Xco likely depends on multiple physical parame- 
ters, such as metallicity, radiation field, gas mass sur 



face density, and density structure (lArimoto et al 



1996 



2006 



Miho s et Pal] H 999: goselli et al] 120021 : iBell et al.l 

Narayanan et al.l 1201 lh ■ Moreover, direct measurement 
of Xco is observationally challenging: For example, 
the assumption of virial equilibrium and the finite spa- 
tial resolution of giant molecular clouds, in even Local 
Group galaxies, may both lead to infl ated values of Xco 
( Blitz et al.ll2007t IBolatto et alll2008f ) ; see iBolatto et all 
(2008) for a more general review of Aco measurements. 
Keeping these complications in mind, our analysis here 
adopts a simple approach: Combinin g the Galact i c mea - 
surement of X C o = 1-8 ± 0.3 from iDame et alj (|2001[ ) 
with the measurements for M31 (Xco = 3.6 ± 0.3) and 
M33 (Xco = 2.6 ± 0.4) from IBolatto et al.l pOOa i. we 
find a mean and range of Xco — (2-7 ± 0.9) x 10 20 cm~ 2 
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(K kms 1 ) . We assume this value to be representative 
of UGC 463, given that the Milky Way, M31 and M33 are 
arguably the only spiral galaxies with well-resolved ob- 
servations of giant molecular clouds or associations from 
which robust measurements of X<jo can be made. 

Figure Q] provides the 24/im image, converted to £h 2 
in units of Mq pc -2 using the hifim/IcoAV calibra- 
tion from Section [2.6.21 and assuming Xco = 2.7 x 10 20 
cm -2 (K kms -1 )" 1 . Figure 0] provides the azimuthally 
averaged surface density profile £h 2 (R) using the 24/im 
surface-brightness profile from Section 12.6.11 Errors in 
Sh 2 are plotted separately for random and systematic 
components; the former includes errors from the /24/im 
calibration and photometry and the inclination, whereas 
the latter includes the /24/im/-/coAU calibration error 
and range in Xqo ■ This estimate of £h 2 (R) agrees 
with the expectation that H2 is concentrated toward the 
"hole" in the H 1 mass surface densit y, also seen in Figure 
Q]and studies of other galaxies (e.g.. IL081 ). 

Out to 15 kpc (4.2 tin), we find M^/M^j = 3.2, a 
value that is reasonable with respect to direct CO and Hi 
studie s in the literature. In particular. lYoung fc Knezekl 
(1989) find a range of 0.2 and 4.0 for, respectively, 
late- and early-type spiral galaxies, comparable to the 
rang e measured by the mo re recent COLD GASS sur- 
vey (fSaintong e et al.l 1201 ID . However, despite having 
a total jj j that is decreased by ~ 10% com- 

pared to the measurement within 15 kpc ([Martinssonl 
|2011|) . UGC 463 is more rich in molecular gas than the 
mean M^/M^ 1 calculated bv lSaintonge et all ()2011[ ) 

for the COLD GASS survey by approximately twice 
the standard deviation. Similarly, by combining our 
24fim imaging, the /24/im /Ico^V and Xg o values de- 
rived herein, and the H 1 data presented by iMartinssonl 
(|2011f ) for 24 galaxies in the DMS, we find a mean of 
•Mjj, /-Wjjj = 0.48, consistent with the iSaintonge et al~l 
(2011) measurement after accounting for the difference in 
their adopted Xco; UGC 463 has the maximum value of 
Xjj 2 /A4jj j and 21 of 24 galaxies have A4jj 2 /A4 H j < 1. 
Therefore, consideration of the molecular mass compo- 
nent in UGC 463 is relatively more important to our dy- 
namical Y» measurements and the baryonic mass budget 
than the majority of galaxies in the DMS. In Section IST61 
we discuss both the total dynamical mass-to-light ratios 
as well as stellar-mass-only measurements to illustrate 
the effects of the gas-mass corrections. 

2.7. Total Gas Content 

Figured] provides a calculation of the total gas content 
of UGC 463, E gas = 1.4(S H i + £h 2 ); the factor of 1.4 
accounts for the helium and metal fraction. For this cal- 
culation, measurements of Sh 2 have been interpolated 
to the radii of the £hi measurements; different inter- 
polation schemes are used in subsequent sections. The 
random and systematic error have been separated as dis- 
cussed above for Eh 2 • Figure Q] demonstrates a strong 
correspondence between \i! K and E gas . Therefore, if n' K 
is a reasonable tracer of the stellar mass, the correspon- 
dence in Figure @] suggests that the radial distribution of 
the stellar mass is roughly equivalent to that of the total 
gas mass. 



3. ON-SKY GEOMETRIC PROJECTION 

Geometric projection parameters are a fundamental 
consideration for the DMS. In particular, accurate in- 
clinations are required to decompose er* into the vertical 
component, o~ z , using the measured axial ratios of the 
SVE (Section [53)) . and to produce the deprojected rota- 
tion curve, V I0 t(R), used in our mass decomposition. Er- 
rors in a z and V T0 t have opposite trends with inclination 
such that intermediate inclinations (25° < i < 35°) are 
preferred and, therefore, used to select optimal galaxies 
for our survey (Papers I and II). 

We present a detailed discussion of our determination 
of the geometric projection of UGC 463 below. We de- 
termine the inclination in Section 13.1) we measure the 
on-sky pointing of each observation in Section 13.2) and, 
in Section [3. 3[ we discuss the two-dimensional maps pre- 
sented in Figure []] created using our final pointing geom- 
etry. 

3.1. Inclination 

We measure inclination using two methods: (1) "kine- 
matic inclinations" (ikm) are determined by modeling 
an observed velocity field by a circularly rotating disk 
and (2) "inverse- Tully-Fisher inclination s" (ztf) are de- 
term ined by inverting the TF relation (|Rix fc Zaritskvl 
1995). Kinematic measurements are used in both incli- 
nation estimates, but in different ways: For itf, one 
measures a fiducial velo city from the rotation curve 
(e.g., Knax or Vfl a t; I VOID in projection and compares 
with the inclination-corrected rotation speed predicted 
by the TF relation for a known absolute magnitude. In 
contradistinction, iki n measurements are independent of 
any distance or photometric measurement, instead de- 
termined by minimizing the difference between measured 
and model isovelocity contours. 

We demonstrated in iPaper iJ that the combination of 
ikin and «tf is ideal for minimizing the errors at low (us- 
ing ztf) and high (using ikin) inclination; the two meth- 
ods produce roughly equivalent errors at i ~ 30°. A com- 
parison of ikin and itf allows for an internal assessment 
of the accuracy and precision of each. We present mea- 
surements of ikin and itf for UGC 463 below, and find 
that the measurements are consistent at 1.2 times the 
combined error, which is satisfactory for our purposes. 
A statistically rigorous combination of the two inclina- 
tion estimates is derived by Andersen & Bershady, in 
prep; however, here we simply produce the error- weighted 
mean value i = 27° ± 2°, which is used in our analysis 
in Section [3.21 and thereafter. One can also estimate in- 
clination via eccentricity measurements of isophotal con- 
tours; however, this method is particularly poor at low 
inclination and for galaxies that have significant outer- 
disk spiral structure, as is true of UGC 463 (Figure [IJ . 
Nevertheless, we calculate a mean isophotal inclination 
of 27° ± 3° by combining Source-Extractor eccentricity 
measurements in the SDSS g, r, and i bands and the 
2MASS J, H, and K bands; this photometric measure- 
ment is easily consistent with our adopted inclination 
based on kinematic measurements. 

3.1.1. Kinematic Inclination 

We use the method described in lAndersen et al.l 
(2008) to measure kinematic inclinations; see also 
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lAndersen fc Bershadvl (|2003f >. The strength of this 
method is in its simultaneous use of the full two- 
dimensional information in our observed velocity fields. 
However, it assumes a single set of geometric projection 
parameters for the entire disk (a "one-zone" model), as- 
sumes that all motion is purely circular rotation in the 
disk plane, and adopts a parameterization for the pro- 
jected rotation curve; therefore, one must justify these 
assumptions. 

Based on edge-on galaxies, literature studies have re- 
peatedly found that warps in non-interacting galaxies 
only influence disk mor phology at lar g e rad ii. Highlight- 
ing two recent studies, Ivan der Kruitl (|2007f) have shown 
that gas disks (as traced by HI) typically begin to warp at 
approxi mately the outer tr uncation radius of the stellar 
disk and iSaha et al.l (|2009f ) have shown that stellar-disk 
warps occur at R > 3/ifl. Thus, we expect no warp- 
ing within the FOV of our IFS of UGC 463, and only a 
marginal warping of our HI data. Indeed, our H I data 
only begin to show a position- angle warp for th e last 
measured radial bin (R = 45": lMartins sonll2011| ). For 
our IFU data, post-analysis of the velocity-field resid- 
uals demonstrates little to no radial dependence of the 
geometric parameters, as determined by translating ve- 
locity residuals with respect to our nominal model (as 
developed in this and subsequent sections) into model- 
parameter residuals. This is done by holding all but one 
parameter fixed and adjusting the free parameter until 
the velocity residual is nearly or identically zero. We find 
no correlation between the parameter residuals and ra- 
dius, with the possible exception of the position angle. 
There is some indication of a positive slope in position 
angle with radius; however, the magnitude of the posi- 
tion angle change is small (less than 5° over the full radial 
range) and the significance of the slope is marginal. This 
means that the use of a radially dependent position angle 
is only marginally justified and, more importantly, in- 
consequential to our measurement of the rotation curve. 
Therefore, a one-zone velocity-field model provides an 
adequate description of our optical kinematic data. 

As briefly noted in Section 13.31 the isovelocity con- 
tours in our velocity fields appear to show slight non- 
circular motions. These motions are most prevalent for 
the Ha data where some coherent structure is seen in 
the velocity-field residual map, particularly along the 
spiral-arm to the south-west of the galaxy center (on the 
approaching side of the velocity field). These coherent 
residuals likely represent streaming motions along this 
spiral arm toward the galaxy center given their spatial 
correlation to the photometric feature and the sign of the 
residual. The magnitude of this streaming is less than 10 
kms -1 along the LOS (less than 25 kms -1 in the disk 
plane) and the covering fraction of all non-circular mo- 
tions is small. Therefore, we expect that our best-fitting 
velocity-field model should suffer only marginally from 
these motions, particularly given the benefits afforded 
the one-zone model in this respect. 

The parameterization of V£?° S (R) = V Ia t(R) sini does 
not adversely bias th e derived geometric param eters: 
lAndersen fc Bershadvl (|2003l ) and lAndersen et al.l (|2008h 
have chosen a hyperbolic tangent (tanh) function — a 
simple two-parameter model that enforces an asymptoti- 
cally flat rotation curve. Although inappropriate for the 



Table 4 

Kinematic Geometry 



Data Set 


*kin 

(deg) 


<t>o 

(deg) 


Vsys 

(kms -1 ) 


Ha 
[Oiii] 
Stars 
Hi 


24.il3:5 

25 ^ +i ' 7 
zo -°-10.8 


68.8+Q3 
68.8+g : * 
68.4l° : ? 
68.8 ± i.5 


4458.6± ;° 
4460.4to.g 
4461.3tof 
4459.5 ± i.5 


Mean 


25.1 ±2.5 


68.8 ±0.3 


4460 ± 1 



rare declining rotation curve in the DMS sample, kine- 
matic inclinations derived for such g alaxies using a m ore 
appropriate parameterization (e.g.. iCourteaul Il997l ) are 
within the formal errors of those derived using a tanh 
model. 

Given our highly sampled and high-quality velocity 
fields of UGC 463, here we use a step function to de- 
fine V^°\R) 1 effectively fitting a set of co-planar "rings" 
with constant rotation speed. We fit up to 13 rings, each 
with a width of 3" (approximately the diameter of a sin- 
gle PPak fiber) such that rotation-speed gradients within 
each ring are small, except for possibly the central ring. 
The final ring includes all data at R > 36" and may be 
omitted for some tracers if no data exist at these radii. 
Despite our use of the term "ring" here, we emphasize 
that this fit is not a typical tilted-ring fit given that we 
are defining only a single set of geometric parameters. 

We measure independent kinematic inclinations for the 
three kinematic data sets provided by the SparsePak Ha 
data and the PPak [O m] and stellar data. All geomet- 
ric and rotation-cur ve parameter s are fi t simultaneously, 
with one exception: IMartinssonl ([201 ID has used recon- 
structed continuum images to determine the morpho- 
logical center of UGC 463 relative to the PPak fibers, 
to which we affix the dynamical center when modeling 
these data. Greater detail regarding our velocity- field 
fitting approach is provided in Appendix [Bj including a 
full description of which measurements are omitted from 
consideration during the fit. However, Appendix [B] is 
primarily focused toward an assessment of the optimal 
data-weighting scheme for modeling the velocity field of 
UGC 463. The rein, we use boot strap simulations (see 
Section 15.6.2 of lPress et al. 2007) to produce inclination 
probability distributions based on four different weight- 
ing schemes. We thereby demonstrate that we obtain the 
most correspondent inclinations among the different data 
sets by adopting weights defined by the derivative of the 
model LOS velocity, Vlos, with respect to the inclina- 
tion, i.e. dV^os/di- These weights approximately follow 
a sin 2 (28) function in azimuth and a direct proportion- 
ality in radius; therefore, data with the most leverage on 
the fitted i nclination (at approximately ±45° from the 
major axis; lAndersen k. Ber shadv 2003) have the high- 
est weight. The best-fitting inclination, position angle, 
and systemic velocity for each tracer are given in Table[4j 
bootstrap simulations are used to calculate the 68% con- 
fiden ce intervals. The r esults provided for our H I data 
from [Martinsson (2011) are based on traditional tilted- 
ring fitting (|Begemanlll989f) . 

The geometric parameters listed in Table [4] for each 
dynamical tracer are in general agreement; the systemic 
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Table 5 

Projected Rotation 
Speed 



Data Set 


(kms- 1 ) 


Ha 


107 ±1 


[Oin] 


108 ±2 


Hi 


105 ±2 


Mean 


107 ±2 



velocities exhibit the most statistically significant differ- 
ences. Such differences are likely due to systematic er- 
rors in the heliocentric velocities of the template stars 
and/or shifts in the pointing center. In any case, these 
shifts are small and irrelevant to our analysis of the 
mass distribution in UGC 463. Using the half width 
of the 68% confidence interval from Table 2] as the error, 
we calculate error-weighted means of ikin = 25? 1 ± 2?5 
and (j>o = 68? 8 ± 0?3. The unweighted mean value 
V S ys = 4460 ± 1 km s _1 has been used in Section 12.11 
to calculate the distance to UGC 463. 

3.1.2. Inverse Tully-Fisher Inclination 



Following the discussion in lPaper III . inverse-TF incli- 
nations are calculated according to 



l T F 



2u r r j dbg 



ci,a - M x 

C2,A 



(3) 



where ci t \ and C2,a are, respectively, the zero-point and 
slope of the TF relation in wavelength band A and 
Mx is the total absolute magnitude. Combining D = 
59.67 ± 4.15 Mpc (the error here is the quadrature sum 
of the random and systematic error from Section 12. ip . 
ttlk = 9.32 ± 0.02 (Section 12.2.3]) . and a /C-correction of 
0.035 mag (|Bershadvll 19951 ). we find M K = -24.59±0.15 
for UGC 463. We use the if-band TF relations derived 
by iVOll to calculate itf based on measurements of the 
projected rotation speed. 

We measure the projected rotation-curve for all gas 
tracers in UGC 463 for use in calculating ?tf; stellar 
measurements are not considered due to significant asym- 
metric drift fSection |4.3|) . Figure [6] presents V^\° 3 for 
the Ha and [O in] data resulting from all four weighting 
schemes implem ented in Append ix [Bj H i measurements 
are directly from iMartinssonl (|2011| ) . It also provides the 

error- weighted mean measurements (V5L_) for data at 
R > 24". No beam-smearing (Section 14. 1[) or pressure 
(Section I4.3.1|) corrections have been applied; these are 
negligible considerations for the measurement of (V^° 3 ). 
The Ha rotation curve exhibits less dependence on the 
applied weighting than does the [O in] rotation curve; 
however, they both compare well with each other and 
with the H I rotation curve, regardless of the weighting 
scheme. The Ha data, in particular, appear to asymptote 
at R > 24"; hence this radial region is chosen for mea- 
suring (V^° 3 ). Table [5] provides (V^° 3 ) for each tracer; 

measurements of {Vj^]. 03 ) from Ha and [O in] are the un- 
weighted mean of the results from all weighting schemes. 
Using all tracers, we find a mean and standard deviation 
of (V r P™ j ) = 107 ±2 kms^ 1 . 



The calculation of itf and its uncertainty relies on 
two additional factors: (1) the choice of the TF rela- 
tion and (2) the estimation of its intrinsic scatter. IVOll 
created multiple samples based on a rotation-curve- and 
asymmetry-based taxonomy, each sample yielding differ- 
ent TF coefficients (ci.a and C2,a) and intrinsic-scatter 
estimates. UGC 463 exhibits Hi and ionized-gas proper- 
ties t hat a re most consistent with the "RC / FD" sample 
from |V0U ( see his Sections 4 and 5 for a detailed defi- 
nition of this sample); UGC 463 fits within the defini- 
tion of the "RC/FD" sample because its rotation curve 
asymptotes to a nearly constant rotation speed and ex- 
hibits neither strong rotation asymmetries (Se ction I4.2p 
nor signs of ongoing interaction. Additionally, IVOll pro- 
duced TF coefficients based on two fiducial velocities, 
one measured at the rotation-curve peak (K nax ) and the 
other where it "flattened" to a constant value (Vfl a t)- 
UGC 463 exhibits a well-defined U nat (Figure [6]) ; how- 
ever, in most cases Vjjat = Vma* for galaxies studied by 
IVOll The smallest obser ved scatter in the if-band TF 
relations derived by IVOll (0.26 mag) was found by ex- 
cluding the outlying NGC 3992 measurements from the 
"RC/FD" sample (leaving measurements for 21 galax- 
ies) and using Vfl a t for the fiducial rotation measurement; 
this TF relation is consistent with having zero intrinsic 
scatter. 

Given the range in c\ t K and ci,k for the "RC/FD" 
sample (with and without NGC 39 92 an d based on ei- 
ther U max or Vfl a t) from Table 4 oflVOl we find 28° < 
*tf < 30°; and we find e(iTF) = 1° and 2° assuming, 
respectively, 0.0 and 0.2 magnitudes of intrinsic TF scat- 
ter, regardless of the assumed c\ } k and ci,k- Taking a 
mean across the four relevant TF relations and assuming 
0.2 magnitudes for the intrinsic TF scatter, we measure 
«tf = 29° ± 2° for UGC 463; this measurement of «tf is 
dominated by systematic error with roughly equal contri- 
butions from the uncertainties in H$ and the TF relation. 
Our conservative approach to measuring itf is justified 
given that, at M K = -24.59 ± 0.15, UG C 46 3 is more 
luminous than any galaxy considered by IVOll . far away 
from the "pivot" point of the fitted F-band TF relations 
(M K ~ -22). 

3.2. Position Angle and Dynamical Center 

The kinematic position angles derived in Section l3.1.1l 
are very consistent among all dynamical tracers; Table 
|4] provides an error- weighted mean value of 0o — 68.8 ± 
0.3, which is constant across the optical disk to good 
approximation (as discussed in Section [3. 1.1[) . As stated 
above, the PPak data affix the dynamical center to the 
morphological center. 

We determine the pointing of each SparsePak IFU ob- 
servation relative to the dynamical center by fitting the 
kinematic geometry (as in Section 13.1. lj) with i and <f>o 
fixed. We simultaneously fit all kinematics measured 
from our IFS, assuming all tracers are in co-planar rota- 
tion. To do so, we apply slight offsets to V sys according 
to the differences found in Table [4] such that all data 
can be forced to have the V sys measured for the stars. 
We also allow for asymmetric drift between the gas and 
stars by simultaneously fitting different rotation curves 
to these components. Finally, we force the [O in] and 
stellar kinematics determined from the SparsePak obser- 
vations to have the same pointing center. During the 
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Figure 6. The projected rotation curves for all UGC 463 gas tracers determined by four different weighting schemes (see Appendix [Fill : 
from light-gray to black, the weights are uniform, error-based, from cos 8, and from dVLOs/di. Left — The projected Ha (filled points in the 
top panel), [Om] (filled points in the bottom panel), and Hi (open points) rotation curves. Measurements of (V^° 3 ) are made at R > 24", 
delineated by the vertical dotted line. The mean (V^" 3 ) for all tracers is plotted as a solid horizontal line. Right — Measured (V^° 3 ) for 
the Ha- (abcissa) and [O lll]-emitting (ordinate) gas for all four weighting schemes; the dotted lines gives the 1:1 relation between the two 

tracers; (V^° 3 ) for Hi is plotted as an open circle on this relation. The black error bars show {V^° 3 ) = 107±2 determined using all three 
gas tracers. 
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with the dither pattern, if allowing for ~ 1" systematic 
errors in the initial pointing. Moreover, reconstructed 
continuum images that use this pointing geometry are in 
good agreement with direct images from SDSS (Figure 
ll Section OS). 



Notes. Columns are: (1) pointing description; (2) nom- 
inal RA dither position; (3) fitted RA position; (4) nomi- 
nal DEC dither position; and (5) fitted DEC position. All 
coordinates are sky-right and relative to the dynamical 
center. 

a The PPak coordinates are taken from Martinsson (2011). 

fitting procedure, we omit velocity measurements based 
on the measurement error and the discrepancy with the 
model, as described in AppendixlBl a nd weight according 
to the velocity errors, as done in lAndersen &; Bershadvl 



Table [5] provides the resulting pointing coordinates rel- 
ative to the dynamical center for each IFS observation; 
as with the geometric quantities in Table El errors are 
68% confidence limits determined using bootstrap simu- 
lations. Table|5]also provides the nominal expectation for 
the pointings based on the dither pattern used during the 
observations. The kinematic fitting results are consistent 



3.3. Two- Dimensional Maps 

Five of the nine images in Figure [T] have used an inter- 
polation algorithm to smooth over the interstitial regions 
of our IFS. The continuum surface-brightness maps of 
our IFS (labeled MjJq, and fJ-Mgi) are determined via a 
calibration to SDSS imaging data. The detailed proce- 
dures used to both perform the surface-brightness cali- 
bration and two-dimensional interpolation are discussed 
in Appendix [X] These interpolated kinematic maps are 
purely for illustration purposes, useful for qualitative as- 
sessments of our registration of the dynamical center and 
a discussion of the two-dimensional kinematic morphol- 
ogy; however, all quantitative analyses herein have been 
performed using the direct fiber measurements, the IFU 
astrometric tables, and our derived pointings. 

The first column of Figure [T] demonstrates the excel- 
lent agreement among the reconstructed continuum im- 
ages and the direct SDSS g-band image. Indeed, the 
central contour of both MjJq, and /J,Mgi directly over- 
lap and are centered on the NED-provided coordinate 
of UGC 463. The detailed spiral structure is apparent 
in, particularly, the //Mgi image due to the small PPak 
fibers. The isovelocity contours of the gas data appear 
to exhibit streaming motions associated with the spiral 
arm toward the south-west of the galaxy center; this is 
less apparent in the stellar data. Additionally, the effect 
of asymmetric drift is seen in the stellar velocity field 
as the "linearization" of the isovelocity contours toward 
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the galaxy center, which is due to a shallower increase 
in the stellar rotation curve, corresponding to the steep 
decrease in the stellar velocity dispersion, toward larger 
radius. We further explore the kinematic axisymmetry 
in Section l4~2l 

4. AZIMUTHALLY AVERAGED KINEMATICS 

Analyses in Sections [5] and [6] assume UGC 463 is ax- 
ially symmetric, considering only the azimuthally av- 
eraged kinematics that we derive in the following sub- 
sections. We apply beam-smearing corrections in Sec- 
tion 14.11 (bea m-smearing corre ctions for the H I data are 
described by iMartfn sson 201 1|) such that kinematic data 
from different instruments can be combined. In Sec- 
tion 14. 2\ we assess the degree of dynamical symmetry 
by comparing approaching- and receding-side kinemat- 
ics. Finding no substantial asymmetries, we discuss the 
azimuthally averaged kinematics in Section 14.31 

4.1. Beam-Smearing Corrections 

Our beam-smearing corrections require a characteriza- 
tion of the beam profile, the convolution of the point- 
spread function and the fiber aperture. Since no sig- 
nificant jitter was detected among or during the indi- 
vidual IFS observations, the effective fiber aperture is 
given by the plate-scale (yielding 2'.'7 for PPak and 4'/ 7 
for SparsePak) . In Appendix El we find that the see- 
ing of the SDSS imaging data — l'/5 in g-band and l'/2 
in r-band — is very close to the effective seeing of our 
SparsePak IFS; iMartinssonl (|2011f) provides a direct see- 
ing measurement of l'/7 for our PPak IFS. Our beam- 
smearing corrections change negligibly over the range of 
measured seeing; therefore, we simply adopt l'/5 seeing 
to calculate all beam-smearing corrections. 

Our approach to beam-smearing corrections (|Paper 11 
depends on comparing our UGC 463 data to mod- 
els of the intrinsic surface-brightness (I), velocity (V), 
and velocity-dispersion (a) distributions. SDSS imag- 
ing data provide the model surface-brightness distri- 
bution; g-band data are used for Mg i-region IFS 
and r-band are used for Ha-region IFS. We assume 
a polyex p arameterization to model t he intrinsic rota- 
tion curve (|Giovanelli fc Havnesl l2002t) . The gas veloc- 
ity dispersion is assumed to be constant with radius 
and isotropic. Only the Ha data are used to describe 
the velocity-dispersion profile; beam-smearing correc- 
tions are marginally different if the [O in] dispersions are 
used instead. For the stars, we adopt SVE axial ratios of 
a = a z /o-R = 0.6 and j3 = ae/a g as determined by the 
epicycle approximation (Section 15.21 equation [5]) . The 
model radial profile for the azimuthally averaged c* (oT) 
combines an exponential function with a cubic Legendre- 
polynomial perturbation at small radius; although some- 
what ad hoc, this form allows for deviations from a nomi- 
nal exponential while enforcing a well-behaved, exponen- 
tial form at large radius. For UGC 463, deviations of aZ 
from an exponential form are small and irrelevant to the 
calculated beam-smearing corrections. 

Beam-smearing corrections are calculated as follows: 
A fit to the uncorrected data is used to generate a seed 
model of the intrinsic galaxy kinematics, which is then 
"observed" by integrating a set of Gaussian line pro- 
files, defined by (I,V,a), discretely sampled over the 



beam profile of each fiber to create a synthetic data set 
(|Westfalll [2009). The velocity and velocity-dispersion 
corrections are the difference between this synthetic 
dataset and the model value at the center of the fiber, 
and they are primarily correlated with the velocity gra- 
dients across the fiber face. The beam-smearing effects 
are largest toward the galaxy center where the rotation 
curve is most steeply rising and the azimuthal coverage 
of each fiber is largest. The trend of the correction is to 
increase the measured rotation speed and decrease the 
measured velocity dispersion. We converge to a set of 
beam-smearing corrections iteratively by updating the 
model of the intrinsic galaxy kinematics, done by fitting 
the corrected observational data, and minimizing the dif- 
ference between the observed and synthetic data sets. 

Monte Carlo simulations demonstrate that the random 
errors in the beam-smearing corrections are < 10%; sys- 
tematic errors, estimated by calculating beam-smearing 
corrections using SVE-shape extrema, are typically much 
smaller. Therefore, we adopt the quadrature sum of a 
10% random error and the estimated systematic error for 
each fiber as the error in the beam-smearing correction; 
the error is always 10% for th e gas dat a. Although lower 
than the upper-limit used in lPaper II this reduction in 
error has little effect on the error budget. 

The correspondence of the uncorrected Ha, [Oili], and 
H I rotation curves in Figure [U despite the factor of > 3 
difference in the beam size among the data sets, suggests 
beam-smearing corrections should be small; this expecta- 
tion is in agreement with our direct beam-smearing cal- 
culations. For the ionized gas data, corrections to Vlos 
are less than 2 kms -1 for 93% of the data, with a maxi- 
mum correction of 14 kms -1 . For the stellar data, Vlos 
corrections are less than 2 km s _1 for 84% and 99% of 
the SparsePak and PPak data, respectively; the maxi- 
mum correction is 7 kms -1 for SparsePak and 8 kms -1 
for PPak. Corrections to er* are less than 5% for 91% and 
99% of the SparsePak and PPak data, respectively; the 
maximum correction is 41% for SparsePak and 29% for 
PPak. Corrections to <7* are typically less than e(<7 bs), 
with the only exceptions occurring near the galaxy cen- 
ter. 

4.2. Axial Symmetry 

Figure [7] presents individual- fiber kinematics after cor- 
recting for instrumental-broadening and beam-smearing, 
with point types indicating the tracer and instrument. 
For measurements located at in-plane azimuths within 

±60° of the major axis, we deproject Vlos to Kot° J i 
rotation velocities are plotted regardless of whether or 
not they were rejected from the velocity-field fitting dis- 
cussed in Section [3] Velocity-dispersion data include all 
measurements made at any azimuth. We find the kine- 
matic measurements from the different IFUs to be very 
well matched. Data are separated according to the ap- 
proaching (negative radii) and receding sides. Figure [7] 
also overlays mean quantities from either side of the mi- 
nor axis. We determine V^" 1 for the ionized gas and 
stars using the velocity-field fitting procedure described 
in Section 13.1.11 with rejection and error-based weight- 
ing. Errors in V^° J are 68% confidence limits calculated 
using bootstrap simulations. The values of cr gas and TTZ 
are error- weighted means. 
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The overlay of the binned data in Figure [7J from its 
two sides show that UGC 463 exhibits little kinematic 
asymmetry, justifying our assumption of axisymmetry in 
the following sections. In detail, the ionized gas rota- 
tion curves exhibit the strongest asymmetry at R < 5". 
UGC 463 is morphologically classified as an SABc galaxy 
(|Paper J ), suggesting that this low- level asymmetry may 
be due to non-circular motions imposed by the presence 
of a weak bar. This kinematic asymmetry may also be 
reflected in the stellar data at marginal significance. The 
velocity-dispersion profiles for both the gas and stars are 
very symmetric at all radii, more so than the rotation 
velocities. 

4.3. Radial Kinematic Profiles 

Figure [8] provides the azimuthally averaged kinematics 
analyzed in Sections [5] and [5J following the same proce- 
dure as described in the previous section but over all az- 
imuth. Stellar kinematics combine both SparsePak and 
PPak observations, and ionized-gas kinematics incorpo- 
rate all tracers from both instruments. Rotation-velocity 
errors (68% confidence limits) are determined using boot- 
strap simulations, not by, e.g., considering the difference 



in rotation speed between the two sides of the rotation 
curve; that is, we assume the disk contains no asymme- 
tries such that any asymmetries manifest themselves as 
an increased error in the measured rotation speed via 
bootstrap simulations. The ionized-gas rotation velocity 
at the largest radius is averaged with the radially over- 
lapping H I measurement to slightly extend the radial 
coverage. The circular-speed curve provided in Figure 
[8] results from applying pressure corrections to the gas 
rotation curve, as described in the next section. 

Measurements of ef^J include only the ionized-gas kine- 
matics, not the colder H I. Although the [O in] veloc- 
ity dispersion is significantly larger than the Ha velocity 
dispersion and decreasing with radius (Figure [7]), the az- 
imuthally averaged cr gas is very nearly the same as the 
Ha velocity dispersion. This is because of the error- 
weighting and the significantly higher quality of the Ha 
velocity dispersions. We find (o- gas ) — 16.6 ± 1.1 kms -1 
when excluding the datum near the galaxy center, which 
is nearly constant as a function of radius. This result 
is comparable to similar measurements made in other 
face-on spiral galaxies by I Andersen et all (|2006fl and, as 
shown by these authors, dominated by turbulence given 



The DiskMass Survey. IV. 



17 







R/h R 
2 





250 


1 

w 


200 


km 


150 


*-> 



u 


100 


> 


50 



1 


50 


w 




(km 




o 
j 




b 


10 



n — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



o V c : Gas 
• V e : Stars 



o cr 



gas 



t^-i^-ns 1 — 



10 20 30 
R (arcsec) 



40 



Figure 8. Kinematic profiles for UGC 463 averaged over all 
azimuth. Top — Circular speed (open circles) and stellar tangen- 
tial speed (filled points), as described in the text. Bottom — Gas 
(open circles) and stellar (filled points) LOS velocity dispersions. 
The dashed gray lines are the best-fitting linear and exponential 
functions for the gas and stellar data, respectively. 



the expected thermal pressure. Physically, the difference 
between the Ha and [O in] dispersions may be related 
to the different energetics involved in generating the two 
lines; for example, [O in] emission may be naturally bi- 
ased toward more turbulent regions of the ISM. From 
Figure UJ we note that the [O in] dispersion is surpris- 
ingly well fit by a radial profile following 2<t7/3 over a 
large radial range; it is of interest to explore the reason 
for this relationship by comparing with other galaxies. 

4.3.1. Circular- Speed Corrections 
iDalcanton k Stilpl (|2010f l derive 



V Z = V 2 



dlncrL, (ilnS 



gas I d\nR 



d\nR 



= V Tot + <Wgas! 



(4) 



where cr gas is assumed to be dominated by turbulence 
and produced by an isotropic gas velocity ellipsoid (cf. 
lAgertz e~a l. 2009). We, thereby, correct our gas rotation 
curve to the circular speed using our measurements of the 
in Figure [8] and £ gas in Figure |U 
One expects cr gas and E gas to decrease with radius such 



that S P > 00 IDalcanton fcStiEl ([2010) propose expo- 
nential functions for use in calculating Sp, which is ap- 
propriate for our S gas measurements; however, we adopt 
the linear function plotted in Figure|8]for cr gas . We calcu- 
late 0.2 < Sp < 3.3 such that the circular-speed correc- 
tions range from 0.4 < (V c — V xo t) < 1-6 kms -1 , always 
below the measurement error in the gas rotation speed. 
If we treat the circular-speed corrections independently 
for [O in] and Ha, we find the Ha-based and [O ni]-based 
circular speeds to be more consistent than if one correc- 
tion is applied to both. However, the combined circular- 
speeds are roughly independent of whether or not the Ha 
and [O in] data are treated separately. 

5. THE DISK 

Using the data described above, we measure physical 
properties of the disk of UGC 463. In summary, we de- 
termine the kinematic scale length, h a , defined as the e- 
folding length of ~&Z (SectionETJ ; we determine the shape 
of the SVE such that we can calculate o z based on our 
measurements of o7 (Section l5.2[) : we determine the total 
dynamical disk mas s surface density, Edyn, using equa- 
tion 9 from lPaper III fSection l5.3p ; we calculate the stellar 
mass surface density, £*, by removing contributions to 
Sdyn from atomic- and molecular-gas (Section 15. 4p ; we 
calculate the stability of the isolated gaseous and stel- 
lar disks, as well as a quantity for the multi-component 
disk (Section 15. 5[) ; and, finally, we measure the dynami- 
cal and stellar mass-to-light ratios in if-band, T^ 1 * (Sec- 
tion ESJ. 

5.1. Kinematic Scale Length, h a 

Measurements of ~oZ are well fit by an exponential pro- 
file; the best-fit exponential has a central dispersion of 
74.7i2.4kms" 1 and h„ = 31'.'3±1'.'8 as plotted in Figure 
[H Given that h a ~ 2.6hp, either T* (or more appropri- 
ately the dynamical disk mass-to-light ratio, Td yn ) or the 
scale height, h z , may be increasing exponentially with ra- 
dius; the implied e-folding length is hph a / ' [he — 2hp) = 
57'.'5. Over the radial range of our data, this sug- 
gests a factor of ~ 2.3 increase in either Td vn (consis- 
tcnt with, e.g. HBell k de Jonell200ll: iZibetti et al.ll2009D 
or h z (compare with edge-on galaxy photometry from, 
e.g.: Ivan der Kruit k . Searlel Il981al: Ide Griis k Peletierl 
1997t lNaravan k Jod|2002t feizvaev k Mitronoval I2002t 
Yoachim k Dalcantonl 120061 : iSaha et all I2009D . The ef- 
fects of radial variations in Td yn and/or h z are further 
discussed in Section [576] 



5.2. Stellar Velocity Ellipsoid Axial Ratios, a and f3 

Calculations of disk mass surface density require mea- 
surements of a z . We obtain u z by correcting ~oZ for 
the shape of the SVE, which can be dire ctly measured 
using stellar and ionized-gas kinematics dShapiro et al.l 
[2003t lNoo7dermeer et all l2008t IWestfall 120091) . We de- 
fine the two axial ratios of the SVE to be a = a z /ap, 
and j3 = o~g/o~p. In the limit where the stellar orbits are 
nearly circular, the epicycle approximation (EA) yields 



fi 1 



Pea 



dlnVe 



+ 1 



(5) 



dlviR 

oc R —1 (Sp = 3) , equati on HI reduces to 
a similar equation used by Swatcrs et al. (2003b) following from 
virial theorem arguments. 
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where Vg is the tangential speed of the stars. Also, as- 
suming that UGC 463 is axially symmetric with an SVE 
that is always aligned with the cylindrical coordinate 
axes, one can approximate 



, AD ,(K 2 -^)sin 2 ^tan^/4g +/32 _ 1 
al 7" V K 



(6) 



from the asymmetric drift (AD) equation 
(|Binnev fc Trem ainc 2008), where we use 







1 



2a 2 



: (1 + /? 2 ) 



(7) 



as defined in IPaper l3 and a and P are constant over 
radial regions where dlna 2 /dlnR « —2R/h a . Using 
these equations, each radially binned measurement of V c , 
Vg, and a 2 provides a direct measurement of Pea and 
a = «ad ) assuming that the derivatives d In Vg / d In R 
and d In al /d In R do not strongly deviate from the fitted- 
model expectations; for all measurements we assume 
al = a7 2 . 

In Figure [9ja) , we calculate Pea for each radial ring 
using the parameterized description of Vg shown in Fig- 
ure [71 Systematic errors of greater than 15% are ex- 
pected due to non-circular stellar orbits (jVandervoortl 
I197& iKuijken fc Tremaindll991tlDehnenlll99~9( ) and non- 
axisy mmetric streaming motions in t he disk near spiral 
arms (|Vorobvov fc Theisl [200l [2008h ; such effects will 
dominate the < 2% random error shown in Figure |9^a). 
Adopting j3 = Pea, Figure [9jb) provides oad given the 
observed values of Sad in Figure [9tc) . Note that the 
calculation of chad near the galaxy center provides an 
imaginary result (a 2 < 0) and is not shown in Figures 
[9ta) or[9]Jb). Excluding this datum, the mean values are 
(a AD ) = 0.46 and (Pea) = 0.81. 

We fit equation [6] to our measurements of <5ad in Fig- 
ure [9jc) under the assumption that the SVE shape is 
constant (such that Sad is linear in radius). We note 
that Sad is inversely proportional to a, and Sad pivots 
about Sad =2 — where <5ad is independent of P — such 
that an increase in 8 decreases the slope of Sad{R)- Er- 
rors in a and 8 are estimated using the quadrature sum 
of two quantities: (1) the standard deviation returned 
by a Monte Carlo sampling of the normal distributions 
i = 27° ± 2° and h a = 31'/3 ± l'/8; and (2) the mean 
of the error determined from a set of 500 bootstrap sim- 
ulations performed for each instance of i and h a . We 
find a = 0.48 ± 0.09 and 8 = 1.04 ± 0.22, shown as the 
dashed line in Figure I9j^l errors are shown as light- and 
dark-gray regions for a and j3, respectively. If we force 
P — Pea when fitting <5ad, we find a = 0.44 ± 0.07 with 
a fit shown by the solid line in Figure |H1 

Figure E[a) shows that (/3 E a) = 0.81 and 8 = 1.04 ± 
0.22 are different at the level of slightly more than the 
random errors. Although interesting, the relevance of 
this difference to our subsequent analysis is primarily 
with respect to our conversion of ~aZ to or (see Sec- 
tion 15.51) and a z . Figure WIA) provides individual mea- 

19 Note that we do not use the average nomenclature for these 
measurements of a and /3, unlike (qao) = 0.46 and (/3ea) = 0.81. 
These fitted parameters describe average properties of the disk but 
they are not strictly averages of multiple measurements. 



surements of the ratio a z /a* = (7 cos 2 i) x l 2 using oad 
and Pea- Adopting i = 27° ± 2°, a = 0.48 ± 0.09, 
and P = 1.04 ± 0.22, we calculate the constant value 
a z /a7 = 0.76 ± 0.09, which is identical to (a z /a7) calcu- 
lated using the individual measurements based on oad 
and /3ea when excluding the datum near the galaxy cen- 
ter. Assuming a = 0.44 ± 0.07 and P = Pea is insignif- 
icantly different from a z j~aZ — 0.76 ± 0.09, as shown by 
the solid line in Figure [9jd). Moreover, fitting the az- 
imuthal variation in a<,/~aZ average d over the full radial 
range , aki n to the analys is done bv lGerssen et al.l (|1997l 
I2000D and IShapiro et all (|2003D . provides P - 0.9, con- 
sistent with the error in our measurement based on 5 ad 
and, therefore, insignificant to the calculation of a z /o7. 

5.3. Dynamical Disk-Mass Surface Density, Sdyn 
We calculate Ed yn using equation 9 from IPaper 11 
which is fundamentally Sdyn = a z /irkGh z . We cal- 
culate a z using ~aZ and our SVE measurements above; 
the oblat eness ratio q — hn/h z , given by equation 1 in 
IPaper III , provides the vertical scale height, h z , at the 
distance D; and the effect of t he assumed ve r tical mass 
distribution is quantified by k. Ivan der Kruitl (|1988l) cal- 
culate the value of k for three vertical mass distribu- 
tions: exponential (k = 1.5), sech (k = 1.71), and sech 2 
(fc = 2). As discussed in IPaper III , our nominal approach 
is to adopt a purely exponential disk as a reasonable ap- 
proximation for the composite (gas+stars) density distri- 
bu tion; thus, of the thre e density distributions discussed 
by Ivan der Kruitl (|1988l ) , we are effectively maximizing 
the measurement of Ed yn and, hence, T^ 1 ^. In IPaper ill , 
we suggested a 14% systematic error in k based on the 
range in k among the exponential, sech, and sech 2 dis- 
tributions. For the case of UGC 463, we increase this 
to k = 1.5 ± 0.3, allowing for a relatively massive gas 
disk (low k) or spherical DM halo (high fc); see further 
discussion in Sections 15.6.11 and [S] 

We use the scale length fitted to \j! k to calculate q 
and h z as provided in Table [2j the band used to define 
liR is insignificant to the result. Figure ITOl plots the re- 
sulting Sdyn calculations along with the measurements of 
Shi & nd Eh 2 from Sections ES] and Section l2.6.3l respec- 
tively. Random (including contributions from V'sys, mx, 
i, Iir, ~a~Z, a, and P) and systematic (including contribu- 
tions from q, fc, and D) errors are calculated separately; 
the former are dominated by the error in the SVE and 
the latter have roughly equal contributions from fc and q 
(~ 20 — 25%). For clarity, Figure [T0l shows neither the 
systematic error in Sdyn nor any errors in Ehi and Sh 2 ; 
the systematic errors in Sd yn are comparable to the ran- 
dom errors. Computing the total gas mass surface den- 
sity, we find a very reasonable fit to Ed yn by scaling S gas 
by a factor of 3, akin to the correlation of E gas with fj,' K 
shown in Figure [4] Efji and Eh 2 have been interpolated 
to the radii of the ~a~Z measurements for this comparison. 

5.4. Stellar Disk-Mass Surface Density, £* 

Calculations of T^ ls ^ require a decomposition of Ed yn 
into stellar and non-stellar components. Five mass com- 
ponents contribute to Sd yn : stars, atomic gas, molecular 
gas, dust, and other non-stellar matter. However, we 
account for only the first three listed; the others are sub- 
sumed in the "stellar" mass surface density, £*, for the 
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Figure 9. SVE results for UGC 463. Dark- and light-gray regions, in any panel, show the error range in /3 = 1.04±0.22 and a = 0.48±0.09, 
respectively, as fitted to the 5ad measurements in panel (c); the dashed line, in any panel, results when these axial ratios are constant 
at all radii. The solid line assumes = /3ea and a = 0.44 ± 0.07 (see text). Panel (a) provides /3eAi calculated assuming the model 
stellar rotation curve from Figure Panel (b) provides a ad calculated assuming /3 = /3ea from panel (a) and measurements of <5ad from 
panel (c). Panel (d) provides measurements of cr z /W^ assuming a = a ad and P = /3ea- Dotted lines in panel (d) give the error range in 
a z /o7= 0.76 ±0.09. 



following reasons. First, our data do not provide dynam- 
ically diffcrentiable measurements of the stellar-remnant 
and other DM (baryonic or otherwise) mass. Second, 
the total dust mass of normal star-forming (or even star- 
bursting) spiral galaxies is expected to be less than 10% 
and more typically ~ 1%) of the total hydrogen mass 
Draine et alj |2007): therefore, it should not be dynam- 
ically important to the mass surface density at scales 
relevant to our analysis. 

Our measurements of E* = Ed yn — 1-4(Shi + ^h 2 ) are 
presented in the bottom panel of Figure rrur As with 
Sdynj we calculate random and systematic errors in E* 
separately; compared with Edyn, E* includes additional 
random error from /24/im and Sjji an d systematic error 
from the i24/xm/-^coAU calibration and Xqo- Given the 
modest errors in /24/im and Sjji, the random error in 
E* is dominated by the errors in o~ z ] and, despite the 
large systematic uncertainties in -?24/jm/-fcoAU and Xqo 
(roughly 30% for both), the systematic errors in £* are 
dominated by the errors in k and q. We plot our S* 
measurements against 2E gas and [i' K (Section l2.2.2j) : the 
scale of the plot is such that all £* measurements falling 
on the fi' K profile have T^g ~ 0.2. We find that E»(i?) 
is roughly consistent with both 2E gas and /i' K at the level 
of the random errors. 

5.5. Global Stability, Q 

We calculate the global stability, Q, of the gaseous 
and stellar disks separately and for the composite disk 



using our mass-surface-density and velocity-dispersion 
measurements These calculations are of interest to our 
dynamical study of UGC 463 given our measurement of 
a systematically low Ed yn with respect to a maximal disk 
(see Section [6]): A dependence between s piral-arm mul- 
tiplicity and disk maxima lity is expected (|Toomrelll981l ; 
lAthanassoula et all 119871) based on stability arguments 
within the context of swing-amp lificat ion theory — se e 
reviews by lAthanassoula! JT984) and iSellwoodT (|2010f ). 
Qualitatively, one expects a higher spiral-arm multiplic- 
ity for lower surface density disks of a fixed rotation 
curve. This expectation is in line with our measurement 
of a submaximal disk for and the three-arm multiplicity 
of UGC 463 (see Figure [TJ; however, it should be noted 
that we find submaximal disks for all 30 ga laxies we have 
studied so far (Paper V; lMartinssonll2011l ) , regardless of 
their spiral structure. Stability issues will be discussed at 
greater length based on our full survey stemming from 
the simple analysis given as an example here for UGC 
463. 

Assuming a razor-thin disk, iToomrel (|1964| ) derived 



Q 



T,gas : 



Qt,* — 



7T G E gas 
K Q~R 

3.36 G E* ' 



(8) 
(9) 



20 Subscripts of Q differentiate between each derivation and 
disk component. 
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Figure 10. Mass surface density measurements for UGC 463. 
The total dynamical (XI,}™; open circles), H I (En,; filled trian- 
gles), H2 (Sjj 2 ; open triangles), and stellar (£»; filled circles) mass 
surface densities are plotted as a function of radius. Random errors 
(black error bars) are plotted for both Edyn and S*; systematic er- 
rors (light-gray error bars) are plotted only for £*. As described 
in the text, dashed lines provide 3E gas and 2E gas in the top and 
bottom panels, respectively. The vertical dotted line marks the 
radius beyond which the Sersic profile describing the intrinsic cen- 
tral light concentration (Section 12.2. 2> contributes less than 3% to 
/j,' K . We overplot fi' K (thick gray line) in the bottom panel with a 
scaling such that S* data falling directly on the surface-brightness 
0.2. 



measurements have TjJ 1 * 



where 



.Vc fVr 



R\R 



dVc 
dR 



(10) 



is the epicyclic frequency and c s is the sound speed in 
the gas; such disks sho uld be stable if Qt > 1. Following 
Uog k, Solomon! (|1984) by treating b oth the gaseous and 
stellar disks as fluids. I Rafikovl (|2001h derived a combined 
stability criterion 
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Qr Qr,* 1 



rui 



;R,gas 



1 



(11) 



where Qr,*/Qt,* = 3.36/7T, QR, gas = QT, gas , u> = 
k\(Xji/ k, k\ = 2it/\ is the wavenumber of the instability, 
and r = c s /ijr; this multi-component, razor-thin disk is 
stable if Qr > 1. The finite thickness of disks systemat- 
ically increases the stability over that estimated via the 
equations above; however, finite-thickness corrections are 
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Figure 11. Disk stability measurements for the isolated gaseous 
(open circles) an d stellar ( filled points) disks following from equa- 
tions derived by Toomrc (1964, Qxi t°p) are plo tted in the top 
panel. The multi-component disk stability from I Rafikovl 120011 . 
Qr, bottom) for wavenumbers of k^ = 1 (black squares) and k^ = 2 
(gray squares) are provided in the bottom panel. Values of QT,gas 
and Qt,*i an d Qr for different k\, are slightly offset in radius 
for clarity. Random errors are plotted in black; systematic errors, 
not provided for Qr when k\ = 2, are plotted in li ght gray. The 
vertical dotted line is the same as plotted in Figure [lOl 



small (|Romeolll992l ). 

Figure HT1 provides our measurements of Q for the disk 
of UGC 463; all measurements demonstrate global sta- 
bility. To calculate k, we use the derivative of equation 
Q]to calculate dV c /dR, adopting the model gas rotation 
curve from Figure [7] and the model cr ga s and Sp functions 
discussed in Section 1431 

Given that turbulence dominates the gas kinematics 
(Section 14. 3[) . we calculate QT,gas by replacing c s with 
our measurements of cr gas from the ionized gas£3 There- 
fore, our measurements of QT,gas may represent upper 
limits; the measured (c g as) = 16.6 ± 1.1 km s _1 for 
UGC 463 is a factor of two or more greater than the 
typical turbulent motions seen in the H 1 or H2 gas 
observed in local spir al galaxies (Hitschfeld et al.l [20091 ; 
iTamburro et all [2009). Figure Ql] shows that Qt. gas in- 
creases monotonically for R > 4 kpc. In the range 
2 kpc < R < 7 kpc (0.6 < R/h R < 2.0), we fuuQ 

21 Reliable measurements of the H I velocity dispersion are dif- 
ficult given both the spatial and spectral beam smearing of our 
UGC 463 data. 

22 Quoted errors in mean quantities here and below are (1) 
the random error in the mean and (2) the quadrature sum of the 
systematic error and the standard deviation in the quantity. 



The DiskMass Survey. IV. 



21 



(QT.gas) = 1-83 ± 0.041^2- Thus, in the absence of the 
stellar disk, the gas disk would be globally stable, but 
only by ~ 1.2 times the error. A factor of two decrease 
in dgas would produce an unstable gas disk, when in the 
absence of the stellar component. 

The calculation of Qt,* is determined directly from 
the data shown in Figures E] and [TU1 where we calculate 
ctr = (1.59±0. 17)57 using our measurements of the SVE 
axial ratios (Section We find (Qi>) = 5.2±2.1+!;£ 
between 2 kpc < R < 7 kpc, decreasing to ~ 0.6QT, g as at 
the radial limit of our data. The isolated stellar disk ap- 
pears to be extremely stable, in stark contrast t o nominal 
expec tations (Q ~ 2) bas ed on empirical (e.g.. IBottemal 
I1993D or theoretical (e.g.. lSerlwood fc Carlberdll984D ar- 
guments. Given the mass of the gaseous disk, a mea- 
surement of the composite stability is physically more 
meaningful. 

The multi-component-disk stability measurement Qr 
asymptotically decreases to Qr ~ 2, more in line with the 
theoretical expectations. In detail, we find (Qr.) = 2.1 ± 
0.4±0.7 at R > hp/ for k\ = 1 indicating a globally stable 
disk£3 measurements assuming k\ — 2, also provided in 
Figure [Til reduce Qr by no more than its error. A factor 
of two reduction in er gas decreases (Qr) at R > Jir by 
4% and 20% for k\ = 1 and 2, respectively, such that the 
composite disk should remain stable in this limit. 



5.6. K-Band Mass-to-Light Ratios, T 



disk 
dyn.K 



id T 



disk 
,K 



We calculate dynamical (T^* K ) and stellar (T*J£) 
mass-to-light ratios using our measurements of Edyn and 
£*, respectively, and the K-b&nd surface-brightness 



if sk = dexp[-0.4(/4 -A\-K- Mq ~ 21.57)] 



1 ^7— I COS I 

robs 



(12) 



in £qPC 2 , where 1^ is the observed emission after 
correcting for Galactic extinction, A l K is t he internal 
dust extinction in mag nit udes (Secti on ]|gg|), K = 0.035 
mag is the /C-correction (|Bershadvl ll995). I K ust /lf s = 
0.02 ± 0.01 is the fraction of dust emission in the K- 
band (Section [2~23]l . i = 27° ± 2° is the inclination 
(Section 153]). and M QiK = 3.30 ± 0.04 is the abso lute 
magnitude of the Sun£3 Equation 10 from IPaper III did 
not include some terms in equation 1121 the most signifi- 
cant of which (0.125 ± 0.019 mag arcsec" 2 ) is the face-on 
correction (cosi); such terms negligibly affect our error 
budget. 

Figure Q2] provides our measurements of Y dl5k - 

and T disk 

A4©/£q; the difference illustrates the effect of the gas- 
mass correction. These data are plotted both as a func- 
tion of radius and (g — i)o color. Given our dynamical 
assumptions (e.g., negligible radial forces), Sdyn rnea- 



dyn,K 

£*//f sk , both in units of 



A total stability calculation following Q gas ^ 
il 11991) is 30% smaller than Q R . 



^R,gas 



Our value and error for Mq x are, respectively, the mean and 
standard deviation of measu r ement s compiled from the fo llowing 
literature sourc e s: IWorthevI 111994 ); IBinnev fc Merrifieldl H1998D; 
IBell fc de Jongl (120011) ; IBell et all 12003) ; and IPortinari et al.l 

pooD . 



surements at small radius may be systematically in er- 
ror, especially in the presence of a pressure-dominated 
bulge (as possible given the central light concentration; 
Section f2.2.2[) . a weak bar ( Section 14. 2|) . or a massive 
halo (Section [6]). The vertical dashed lined in Figure IT21 
marks where our model of the central light concentration 
contributes less than 3% to the total light; Sdyn mea- 
surements beyond this radius (i.e., R > 2 kpc) should 
not be strongly affected by either the "bulge" or the dy- 
namical assumptions. The two measurements of T^^ K 

and T^ 1 ^ within this radius are not plotted as a func- 
tion of (g — i)o- Measurements of (g — i)q are interpolated 
from our SDSS photometry (Section l2.2.1[) after account- 
ing for internal reddening (Section I2.2.4P ; instrumental- 
smoothing- and /C-corrections have not been applied to 
the g- or i-band data. We compare our dynamical mea- 
surements with SPS model calculations in Figure [12] by 
plotting 

(13) 



T^ = dexp[oj C + 6j S -(s-*)o], 



an and bx are 
-0.211 and b K = 



taken from 
0.137, here- 
1.379 and 



where the coefficients 
IBell et all (12001 g K = 
after IB03I ) and Izibetti et al.l (|2009l a K 
bK = 0.604, hereafter IZ09f ) ; these models roughly repre- 
sent the extrema of si milar SPS modeling done b y, e.g., 
IBell fc de Jongl (|200l and IPortinari et al. (2004). Con- 
sistent with our photometry in Section |2T2"1 IB03I and lZ09l 
assume i-T-band measurements are in Vega magnitudes 
and g- and ?-band measurements are in AB magnitudes. 
We discuss our results below after generating approxi- 
mate probability distributions for T^ is f(R) and (T^ 1 ^). 

5.6.1. Probability Distribution 

We create an approximate probability distribution 
(combining random and systematic components) for our 
measurements of T^(R) and (T^) by Monte Carlo 
(MC) sampling individual probability distributions for 
each component of our T^ 1 ^ calculation. In contrast to 
the random error contributors, systematic errors — par- 
ticularly in q, fc, and Xco — may not be normally dis- 
tributed, providing the primary motivation for this test. 
Our simulation is limited by the exclusion of any param- 
eter covariance, such as might be manifest in a refitting 
of a and j3 after adjusting i; however, we do not expect 
parameter covariance to dramatically change the funda- 
mental conclusions drawn from this MC simulation. 

Random-error contributors — all quantities with mea- 
surement errors contributing to the random error in T dl * 
- are assigned Gaussian probability distributions ac- 
cording to their derived e. We also assign a Gaussian 
distribution for D (combining the random and systematic 
error in quadrature) and l24fim/ IcoAV (using the ±0.11 
dex error in equation [TJ). Information on the probability 
distribution for Xqq , particularly for spiral galaxies like 
UGC 463, is limited; therefore, we simply assume a uni- 
form distribution with the range X C o = (2.7±0.9) x 10 20 
cm~ 2 (K kms -1 ) -1 (Section 12. 6. 3[) . We h ave derive d an 
empirical probability distribution for q in IPaper III ( see 
Figure 1, and references, therein). For our MC simula- 
tion, we smooth the growth curve of this empirical distri- 
bution by a low-order polynomial to avoid the discrete- 
measurement quantization noise. 
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Figure 12. Measurements of T^ 1 * (filled points) and K (open circles) as a function of radius (left) and (g — i)o color (right); the 

^~dyn K meas urements are slightly offset in radius from the T^ 1 * measurements for clarity. Random errors are plotted in black; systematic 
errors, only plotted for TjJ'|^ in the left panel, are light gray. The vertical dotted line delineates the radius at which the central light 
concentration contributes ~ 3% to the total flux (R = 2 kpc); data within this radius are not shown in the right panel. The solid dark-gray 
line provides the mean Tf g at R > 2 kpc. The (g - i) colors are used to predict Y^ 1 ^" based on the work of B03 (dotted line) and Z09 
(dashed line); random errors due to the photometry are shown by the surrounding gray regions in the left panel. 



A robust probability distribution for k is elusive, lack- 
ing an empirical measurement. For UGC 463, the mea- 
surements of a relatively massive gas disk (Section 12. 7\i 
and DM halo (Section [S]) are particularly relevant to the 
value of A;. In an extreme scenario, the stellar disk is ver- 
tically exponential and the gas disk is razor thin, yielding 
an effective 30% decrease in k for our measured values 
of E gas and £». Despite the resulting increase in Sdyn, 
we would still infer a massive DM halo that increases 
k by 20 — 30% (Section 16.3. 2p and, therefore, roughly 
offsets the effect of the massive gas disk. In view of 
further complications introduced by finite-thickness gas 
disks, multi-component stellar disks, and triaxial halos, 
we have decided to take a simple approach and assume 
k is quantized and equally distributed among the expo- 
nential {k = 1.5), sech (k = 1.71), and sech 2 (k = 2.0) 
cases. 

The results of our MC simulation are shown in Figure 
[13] based on 10 6 recalculations of T^ 1 * for each measure- 
ment of ~oZ at R > 2 kpc. Each recalculation is binned 
in the two-dimensional (R, YjJ 1 ^) plane; R is binned in 
physical units, incorporating the MC sampling of D. For 
each radial bin, we create a growth curve for T^ 1 * such 
that higher intensity (darker) cells in the left panel of 
Figure[13]represent more probable measurement of Y^ 1 ^. 
We overplot 68%, 95%, and 99% confidence contours for 
T^ 1S ^(R), as well as a contour following the median value. 

We also overplot the nominal measurements of T^ ls ^ from 
Figure[12]for reference, again differentiating between ran- 
dom and systematic error. Figure [13] also provides the 
growth curve of (Y^ 1 *), calculated for each of the 10 6 

recalculations of T^ 1S ^(R). 

5.6.2. Discussion 
Based on our nominal calculation (Figure [T2l) . we find 
( T dytx) = °- 34 ± °- 09 ± °- 15 and ( T *]k) = °- 22 ± 



0.09±[J^ at R > 2 kpc. in units ofM /£g. In contrast. 
El and [H predict (T^) = 3.6(T^f) and (T^) = 
0.56(T^ 1 ^), respectively. From our derived probability 
distribution we find a median of (YjJ 1 ^) = O.Ht'om sucn 
that the IZ09I pr edicti on is within our 68% confidence in- 
terval; however, IB03I predict a measurement that occurs 
for less than 1 in 10 5 recalculations of (Y^ 1 *). The me- 
dian value of our probability distribution is below our 
nominal measurement because the median value for k is 
above our nominal value of k = 1.5. Within the pa- 
rameter space probed by our MC simulation, we find a 
maximum measurement of (Y^ 1 *) = 0.82. 

The prescriptions adopted for the many ingredients 
of SPS modeling — such as the initial mass function 
(IMF), star- formation and chemical-enrichment history, 
dust content, and the treatment of specific phases of 
stellar evolution — play an important role in s etting 
and its tren d with color (jPortinari et all 120041 : 
IConrov et all [20091 120101: IConrov k Girnn] [2010h . In- 
deed, in their discussion of the dif ferences b etween their 
Y* -color relations and those from lB03l IZ09I isolate their 
treatment of the star-formation histories and thermally 
pulsating asymptotic giant branch (TP-AGB) phases of 
stellar evolution from intermediate-age populations as 
the primary culprits. The latter particularly effects dif- 
ferences in SPS predictions of Y, in the NIR bands 
(|Marastonll2005D . As ide from the factor of — 7 differ- 
ence in the mean iBOl andUl predictions for UGC 463, 
it is also important to keep in mind that there is an addi- 
tional factor of 2 or more internal variation in the Tf P x 
calibration associated with each study as determined by 
their search of the SPS modeling parameter space. It is 
encouraging that t he ad vancement in SPS modeling, as 
represented by the Z09 study, are more consistent with 
our dynamical measurements. However, Figure [12] shows 
that neither the zeropoint nor the trend of Y* with (g — i) 
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Figure 13. Result of the Monte Carlo simulation outlined in Section 15.6.11 The left panel shows the gray-scale image of the folded 
growth curve in T^ 1 * at a given radius. The median T^'*(i?) lies at the peak intensity (black) as traced by the dashed white line; black 
contours mark the 68%, 95%, and 99% confidence intervals. The nominal measurements of T^ 1 *(i?) from Figure [T51 are overplotted (black 
points) with random and systematic errors in black and gray, respectively; data at R < 2 kpc are omitted. The growth curve of (T^ 1 *) is 
plotted in the right panel. Gray lines mark (Tf B £) = 0.22 and 0.476; less than 1% and 0.1% of the simulations produce (T^) > 0.476 
and (T^*) > 0.59, respectively. 



color from these models are a good match to our mea- 
surements; therefore, we cannot conclude that this spe- 
cific SPS treatment — either for the TP-AGB phase or 
star-formation histories — is correct or applicable to our 
entire sample. An absolute calibration of T* using a 
multi-color approach for the full DMS Phase-B sample 
will be presented in forthcoming papers. 

As expected from the discussion in Section 15. 1[ Fig- 
ure Q2] shows TJjy^ K and T^J 1 ^ generally increase with 
radius. Considering only the data at R — 5.6 and 10.8 
kpc, Tf^ K and Tf$ increase by a factor of 2.2±1.0t^ 

and 3.4±2.33jJ, respectively. Figure fT3l shows that the 
radial trend holds for the median of the probability dis- 
tribution in T^ ls x(R). Variation in T;^ is expected in 

galaxy disks given the observed arm/inter-arm and radi- 
ally averaged color gradients. Yet, in terms of the latter, 
Figure QU shows that the variation in Yf$ for UGC 463 

is consistent with neither of the plotted Tf p ^ predictions. 
Given both the errors in our measurement and the errors 
in the calibration, it is difficult for us to conclude 

that our measurements are inconsistent with a radially 
invariant T^ 1 ^. Indeed, a radially independent T^ 1 ^ is 
consistent with our 68% confidence limits derived in Fig- 
ure Q22 

Moreover, the more shallow decline of <j\ with respect 
to the surface-brightness profile, which is the primary 
driver for our measurement of a radially varying T^ 1 ^, 
may also be interpreted as a flaring of the stellar disk 
(Section l5.ip . A flared stellar disk has been measured for 
the Galaxy (iLopez-Corredoira et al.1 120021 : iRobin et al.1 
120031 : [Momany et al.ll2006| ) and one might expect disks to 
be flared due to, e.g., interactions with dark and/or lu- 
minous satellite s (Hava shi fc Chib a 2006; Dubinski et al.l 
2008tlRead et al.ll2008t iKazantzidis et alJl2009Q . Indeed. 



a nearly constant a z at large radii in M 83 and M 94 
provide evidence for such interactions. However, the 
onset radius for disk flaring is expected to be beyond 
the region relevant to our dynamical measurements for 
UGC 463. Robust photometric evidence for stellar-disk 
flares in edge-on galaxies remains elusive: The empiri- 
cal and theoretical foundation for the vertical structure 
of galaxy dis ks with radially independent sc al e heigh ts 
developed by Ivan der Kruit fc Searld (11981allbl Il982allbh 
remains the current paradigm due to repeated confirma- 
tions of little to no variation in scale height measure- 
ments from surface photometry, particularly for galax- 
ies of similar Hubble type to UGC 463 (SABc) (e.g . , 
Ide Griis fc Peletierl 119971: iBizvaev fc Mitronoval 2002) 



Herrmann et alJ ( 2009( 1 suggest their measurements of 



However, claims of factors of 2 or more increase in scale 
height within the optic al extent of some stellar disks exist 
in th e literature (e.g.. iNaravan fc Joel 120021 : iSaha et al.1 
2009). For our UGC 463 data, the measured increase in 
T^ 1 ^ with radius is consistent with these claims; how- 
ever, we cannot claim a stellar-disk flare exists in UGC 
463 based solely on our data. Therefore, barring more de- 
tailed information on the disk structure of UGC 463, we 
simplify our mass decomposition in Section |5] by largely 
focusing on results that assume T^ 1 * and h z are constant 
for the entire disk. 

Finally, we note that our data provide a few, limited 
assessments of the presence of a "DM disk" in UGC 463. 
Such a structure has be en predicted by recent simulations 
(e.g.. iRead et al1l2009fl and modeling of the Galaxy by 
iKalberlal (|20'03T Wcf.. iMoni Bidin et al.ll2010D . DM disks 
are expected to be more exten ded both r adiall y and ver- 
tically than stellar thin disks; IKalberlal (120031 ) fit a DM 
disk that has a scale length and scale height that are, 
respectively, three and 10 times larger than for the stel- 
lar thin disk. Our stellar kinematic data are expected to 
trace the thin disk mass distribution only such that we 
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can place an upper limit on any DM distributed iden- 
tically to this structure as follows: Assuming the Tf p ^ 

prediction from IZ09I is exactly correct, our Y dl * mea- 
surements suggest a thin DM disk that has ~ 80% of 
the stellar mass surface density (or ~ 30% of the bary- 
onic mass surface density). One can increase the mass 
surface density of a DM disk in UGC 463 by proportion- 
ally increasing its scale height with relatively moderate 
effects on our calculation of Edyn- Assuming no influ- 
ence on Edyn, the scale height of the DM disk would 
need to be roughly the same a s the thin-disk sc ale length 
to reach the mass ratio that [Kalbcrla ( 2003) measure 
for the Galaxy. Although one may accommodate such 
a disk within the current understanding of DM disks, 
such a structure is mostly conjectural with respect to 
our data. Given the uncertainty in both the expected 
vertical distribution of a DM disk and the SPS model- 
ing results for T^ 1 ^, our mass decomposition assumes no 
DM disk exists in UGC 463, which is consistent with our 
observations. 

6. MASS BUDGET 

We produce a detailed mass budget of UGC 463 
via a traditional rotation-curve mass decomposition 
(e.g..lCarignan fe Freemanll 19851 : Ivan Albada et al.lll985t 
iBegeman et al.lll991f ). We assume 

K 2 = £^ 2 > (14) 

3 

where Vj is the circular-speed of a test particle associated 
with potential-density pair $j and pj for each mass com- 
ponent j; all potential-density pairs are considered inde- 
pendent and separabl e, neglecting any covariance amon g 
the j components (cf . . lvAS86t I Amorisco fe Bertidl20"Ioh . 
We calculate Vj for each baryonic mass component based 
on our mass-surface-density measurements using rotmod, 
a program within the GIPSlffi software package that cal- 
culates Vj for oblate and s pherical density distributions 
(following ICasertariolll983l ) . 

Studies of our Galaxy suggest the potential-density 
structure of UGC 463 may be very complex. For simplic- 
ity, we assume here that the total gravitational potential 
is composed of four unique, axisymmetric density distri- 
butions, yielding the following circular speeds: (1) Vhaio 
for the spherical halo, (2) V* bulge for the stellar bulge, (3) 
ydiak f or ^ e stellar disk, and (4) V^ k for the gaseous 
disk. We attribute the central light concentration to a 
bulge (cf. Dutton 2009, who propose such surface den- 
sity peaks may be attributed to the disk); however, our 
use of the term "bulge" does not distinguish between a 
(disk-like) pseudo-bulge ([Kormendv fc Kennicutti 120041) 
or (spherical) classical bulge. We neglect any signifi- 
cant contribution from, e.g., an inner halo, thick stellar 
disk, or flattened dark-matter component for two rea- 
sons: (i) Our spectroscopic and imaging data present 
no evidence for significant contributions of such compo- 
nents to either the measured kinematics or the stellar 
light profile; and, therefore, (ii) if present, the influence 
of such mass components is negligible with respect to our 

25 Groningen Imag e Processing System; 

http : //www . astro . rug . nl/ - gipsy/ 



Table 7 

UGC 463 Enclosed Mass at 15 kpc (i.2h R ) 







Mass Fraction 




Mass 


Baryonic 


Total 


Component 


(10 10 ^1 Q ) 


(%) 


(%) 


Stellar Bulge 


0.22 ±0.09 ±0.16 


5.2 


1.2 


Stellar Disk 


2.6 ± i-i±i;l 


62 


15 


Total Stars 


2.8 ± 1.2± 2 -° 


67 


16 


Atomic Hydrogen 


0.24 ±0.02 


5.7 


1.4 


Molecular Hydrogen 


0.76 ±0.041^* 


18 


4.3 


Total Gas 


1.40 ± 0.071^43 


33 


7.9 


Baryonic Matter 


4.2 ±l.l + i9 




24 


Dark Matter 


10 =+1.9 +3.2 
10 - J -2.4 -4.2 




76 


Total Mass 


-| rt iy-|-2.2 +3.8 
1 ' ■ ' -2.7 -4.6 







analysis of the gravitational potential in UGC 463. Fur- 
thermore, we assume the halo is dominated by DM (i.e., 
Uhaio ~ Vdm) and the vertical distributions of the atomic 
and molecular gas are identical (V^ s sk is determined by 
E gas directly). We define the baryonic circular speed, 
vl = (U* bulgo ) 2 + (U disk ) 2 + (U g d i sk ) 2 , such that we can 
isolate the DM mass contribution via V-£ M = U c 2 — U b 2 . 
The fundamental advantage of our rotation-curve mass 
decomposition over previous studies is that Vb is uniquely 
defined by our observations. 

In summary, we calculate the circular speeds of all 
baryonic mass components in Section 16.11 we mea- 
sure Vdm(-R) an d the DM-halo volume-density profile, 
Pdm(-R), in Section [6.21 and we discuss the relative con- 
tributions of the baryonic matter and DM components to 
the total mass budget of UGC 463 in Section 16.31 finding 
that DM dominates at R > fiR. 

6.1. Baryonic Mass and Circular Speed 

Follow ing th e bulge-disk decomposition discussed in 
Section I2.2.2[ we calculate the mass-surface-density 
distribution for each stellar component assuming 
T* tK (R) = (T d jf } = 0.22 ± Q.09±g:i|. We use a single 
Y*^ for the disk and bulge components, which is rea- 
sonable considering the marginal change in T^ 1 ^ between 
the disk- and bulge-dominated regions in Figure [12] The 
results are shown in Figure [H] along with our measure- 
ments of E gas . We calculate masses enclosed within 
R = 15 kpc (R = 4.2/ifl) for all baryonic components 
in Table [7] by integrating each mass-surface-density pro- 
file; these results are further discussed in Section [5.3.31 
We note that the instrumental-smoothing correction to 
p! K (Section l2.2.2[) amounts to a marginal 3% increase of 
the stellar-bulge mass. 

Circular-speed calculations for each baryonic compo- 
nent use the mass-surface-density profiles in Figure Q3] 
and an assigned three-dimensional, axisymmetric density 
distribution. We assume any truncation of each mass el- 
ement occurs well beyond our last dynamical measure- 
ment. We assume E gas is distributed in a razor-thin disk; 
and, consistent with our previous assumptions, we adopt 
an exponential vertical distribution for the stellar disk 
with a constant scale height of h z = 0.44 kpc. For the 
stellar bulge (or central mass concentration), our nomi- 
nal approach is to assume a spherical distribution. 
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Figure 14. Mass-surface-density (top) and circular-speed (bot- 
tom) distributions for the stellar bulge (dotted line), stellar disk 
(dot-dashed line), gaseous disk (dashed line), and all baryonic mat- 
ter (solid black line) in UGC 463 assuming a constant T, k(R)- 
Our measurements of S^yn at R > 2 kpc are overplotted as open 
symbols in the top panel for reference, which are directly traced 
during our calculations assuming a variable T« k(R). Dark- and 
light-gray regions illustrate, respectively, the random and system- 
atic error in Sb and Vt>- The solid gray lines assume T tj x = TT^ 1 ^ 
as predicted by |B03t Measurements of the DM-only circular speed 
(Vdm; black triangles) are calculated by subtracting Vb from V c 
(open triangles) in quadrature. Random and systematic errors in 
Vdm are shown as black and gray error bars, respectively. 

Due to our imposition of a constant T*k for the 
mass decomposition illustrated in Figure [T4J the sur- 
face density of all baryonic mass components, Sb, does 
not exactly follow our Sd yn measurements. We have 
also performed a more direct mass decomposition that 
adopts Eh = Sdyn by applying a smooth interpolation for 
T*,k{R) = Tf$(R) at 2 < R < 11 kpc; we extrapolate 
by simply extending the Y^ 1 * measurements at the two 
endpoints of this range to small and large radius. This 
approach does not produce statistically different values 
for the quantities discussed below. For completeness, Ta- 
ble [8] provides the implied DM properties for both a con- 
stant and variable ¥* k(R); however, we simplify the 
discussion below by primarily focusing on the results ob- 
tained by assuming a constant Y*j<(R)- This approach 
is justified by the consistency of our measurements with 
a radially invariant T^ 1 ^ as discussed in Section 15.6.21 



In addition to the constant Y* t K based on our dynam- 
ical measurements, Figure [14] also provides the calcu- 
lation of Sb and Vb assuming Y*^ = Tf p ^ = 0.79, 
based on the prediction of B03; adopting the IZ09I pre- 
diction produces Vb within the systematic errors of our 
dynamical measurements. The stellar and baryonic mass 
are, thereby, increased by a factor of 3.6 and 2.7, re- 
spectively, resulting in a maximal disk (i.e., Vb ~ V c ). 
Therefore, we can account for the measured rotation ve- 
locity within R < 3.2/i/j by simply scaling fi' K by a con- 
stant Y*^. This result is not u nexpe cted given that 
the "diet Salpeter" I MF, used by IB03I . was chosen by 
iBell fc de Jongl (|200Tl ) to accommodate the "maximum- 
disk" rotation-c urve mass decompositions produced by 
iVerheiienl (|1997| ). Moreover, this is consistent with the 
expectation from most rotation-curve mass decomposi- 
tions in the literature p erformed wit hin similar radial 
regimes (see discussion bv lSancisill2004h . in cluding the re- 
cent s tudy of dwarf galaxies performed by iSwaters et al.1 
(|2011l) . However, we have shown in Section [5.61 that the 
assumption Y*^ = Yf^f = 0.79 is rejected by our Y^ 1 ^ 
measurements at 3> 99% confidence. That is, while it 
is possible to account for the rotation curve of UGC 463 
within R < 3.2/i# without invoking DM, DM is effec- 
tively required by our treatment of the observed stellar 
kinematics and, in fact, dominates the mass budget (Sec- 
tion ^. 3. 3j) . Moreover, the implied DM-mass distribution 
must be substantially less oblate than the baryonic-mass 
(disk) distribution traced by our stellar kinematics to si- 
multaneously explain our measurements of V c and <t7, as 
discussed in Section ["5. 6.21 

We note here that the DM properties we infer below are 
robust against many of the assumptions made above con- 
cerning the detailed baryonic mass decomposition. First, 
our characterization of the intrinsic central light concen- 
tration is a marginal consideration for calculating Vb- 
For the innermost measurement of V c (R = 0.7 kpc), Vb 
is 10% (9.1 kms -1 ) higher for our nominal measurement 
than if we were to adopt the stellar-disk density distri- 
bution at all radii; the random errors are 15%. This 
difference increases to 25% (23.3 kms -1 ) if we also omit 
the instrumental-smoothing correction to /i' K ; however, 
the difference at the radius of the second measurement 
of V c (at R = 1.5 kpc) in this case is only 2% (3.5 
kms -1 ). Second, we find statistically equivalent mea- 
surements of the DM rotation speed, Vdm , when simply 
calculating Vb directly from Sb = Sd yn and assuming 
the nominal oblateness of the disk. The detailed baryonic 
mass decomposition presented here is meant to couch our 
rotation-curve mass decomposition within the traditional 
construct found in the literature and to provide insight 
into the baryonic mass budget of UGC 463. 



6.2. Dark- Matter- Halo Circular Speed, Mass, and 
Volume Density 

Figure[T4lprovides measurements of V£ M — V^—V^ as- 
suming a constant Y*^ (R); the enclosed mass of the pre- 
sumed spherical halo is M D a ^ = (8.9±{i tH) x 10 10 M Q 
within the radial range of our measurements (i? < 11.4 
kpc). Measurements of Vdm also provide the spherical- 
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Figure 15. The DM-only circular speed (Vdm; top) and volume- 
density profile (pdm; bottom) for UGC 463 assuming a constant 
T t k{R)'i random and systematic errors are shown in black and 
gray, respectively. Upper limit arrows are used for data with ran- 
dom errors of e(pdm) > PDM- The best-fitting NFW (solid lines) 
and pseudo-isothermal (dashed lines) parameterizations of the DM- 
halo are shown against the data. For reference, dotted lines in the 
bottom panel provide simple power-law density distributions fol- 
lowing pdm oc R~ a for a = 1 and 2; zero-points are set to approx- 
imately match the trend at small and large radius for a = 1 and 
a = 2, respectively. 



halo volume-density profile via 



Pdm 



V 2 

'DM 

4irGR 2 



d\nR 



(15) 



(|de Blok et all 120011: ISwaters et all I2003al) as shown in 
Figure [TBI where we calculate dV£ M /dR = dV 2 /dR — 
dV^/dR based on the model values for V c and V\>. The 
calculation of dV c / dR is the same as used in Section 15.51 
to calculate n, and the calculation of dVb/dR is done via 
finite-differencing of the V\> data in Figure [L4j Measure- 
ments of pdm are upper limits at R > 4 kpc (R > l.lha) 
because the random error, c(/5dm), becomes larger than 
the measured value. 

We model t he Vdm data in Figure [T5l using a spherical 
"NFW" halo dNavarro et al.llilM fl997j) and a pseudo- 
isothermal sphe re (Ostrike r fc Caldw ell 1979; Sc hmidt! 
119851 : iKent 1986|). The former is parameterized by the to- 
tal halo mass, A^ D ^, and its concentration c = R2oo/R s , 
where -R200 is the radius at which the halo reaches 200 
times the critical density (p C rit = 3Hq/8ttG) and R s is 
the characteristic scale of the density profile; as above, we 
adopt Ho = 73 kms -1 Mpc -1 . The pseudo-isothermal 
sphere is parameterized by its central density, p$, and 
core radius, R c . Best-fitting DM halo parameters are 
determined by ^-minimization, where \ 2 lS defined us- 
ing only the random errors. Bootstrap simulations are 
used to calculate the random error in each the parameter. 
Systematic errors are based on 100 simulations of refit- 
ting the DM halo after MC sampling of the systematic 



errors in Vdm, again minimizing a \ 2 incorporating only 
the random errors. The best-fitting DM halo models are 
plotted against our measurements in Figure 1151 we pro- 
vide the best-fitting parameters for both a constant and 
variable T^ K (R) in Tabled] Statistically, the NFW and 
pseudo-isothermal halos are equally suitable descriptions 
of the DM halo of UGC 463, and there is only a marginal 
difference between the results when using a constant or 
variable T*,k{R)- While the reduced % 2 is lower when 
adopting a constant Y*jf (i?), the systematic deviations 
of the data about the model is such that this difference 
is not statistically meaningful. For comparison, Figure 
1151 also shows simple power- law density profiles following 
Pdm oc R~ a for a = 1 and 2; zero-points have been set 
by-eye in each case. 

We measure the concentration of the NFW halo to be 
c = 19^ i n reasonable agreement with the expectation 
for low-redshift DM halos of similar mass (jNavarro et al.1 
119971: iBullock et all 120011) . Thus, consistent with our 
measurement of a relatively low-mass disk, it appears 
that the baryons have had little effect on the structure 
of the DM halo. However, it is difficult to assess the 
physical meaning of this concentration index given the 
simplicity of our halo fitting; specifically, we do not in- 
clude adiabatic contraction, which should be a relatively 
small effect given the expected gravitational influence of 
the baryons. It is interesting that the innermost Pdm 
datum indicates a steeper slope than provided by either 
of the DM-halo parameterizations, demonstrating a ~ 2 
as opposed to unity. However, the error in this measure- 
ment is large and is highly subject to our treatment of 
the central mass concentration. The full DMS sample 
will provide stronger statistical constraints on the inner 
halo profile. 

6.3. The Dominant Gravitational Influence of Dark 

Matter 

As discussed in Section [TJ assessments of the relative 
contribution of dark and baryonic matter to disk-galaxy 
mass budgets have been limit ed by the disk-halo degen- 
eracy (jvan Albada et aflll985[ ) . We have measured these 
contributions directly for UGC 463. Although allowing 
for a novel investigation of the DM mass distribution, 
it is useful to cast our analysis also in terms of a tradi- 
tional approach such that we can compare with previous 
studies. 

In particular, a common praxis in rotation-curve mass 
deco mpositi on is the so-called "maximum-disk" hypoth- 
esis (|vAS86l ) , producing an extremum of the mass bud- 
get. Application of the "maximum-disk" hypothesis 
amounts to maximizing a radially independent mass-to- 
light ratio while adjusting any DM-halo parameters to 
fit the observed circular speed; however, exact imple- 
mentations have varied. All baryonic components have 
been approximated by a single exponential disk (as in 
Ivan Albada et al.lll985|). val ues of T* have been distinct 
(as in. e.g.. lKentlll986lll987l ) or identical (as done above) 
for the bulge and disk components; and the gas disk has 
been subsumed into or isolated from the stellar compo- 
nents). Thus, one should keep in mind that, although 
termed the "maximum- disW hypothesis, the direct asso- 
ciation of this hypothesis with the stellar disk, in par- 
ticular, can be tenuous in its practical implementation. 
Regardless, all implementations minimize the DM con- 
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Table 8 

UGC 463 DM-Halo Properties 



Parameter 


Constant Y* x 


Variable k 


Pscudo- Isothermal: log po 


8.85 ±0.21 ±0.13 


8.83 ±0.16 ±0.09 


Rc 


1.06 ±0.21 ±0.01 


1.16 ±0.21 ±0.04 


xl 


0.64 


1.20 


NFWHalo: logXg^° 


11.51 ±0.09 ±0.15 


11.60 ±0.09 ±0.07 


c 


19.2 ± 2.7 ± 1.6 


18.9 ±2.7 ±1.2 




0.61 


1.14 


J=b(-Re) 

•7^(2.27^) 

J r f sk (2.2/i i{ ) 


n H1 +(J.4W +0.24 
u - ol -0.15 -0.35 
n ci +0.07 +0.12 
u -° -0.09 -0.18 
n AR+ - 09 + - 15 
u -™-0.12 -0.22 


n 7 o+().43 +0.17 
u -'°-0.13 -0.22 
n co+0.05 +0.09 
u - do -0.07 -0.12 
n o fi +0.08 +0.11 
U - JO -0.10 -0.15 



Notes. Units of p , R c , and -Mdm° are -M.® kpc 3 , kpc, and Mq, 
respectively. 



tribution to the mass budget at small radius. Galaxies 
adhering to the "maximum-disk" hypothesis are often 
said to have maximal disks; however, this definition re- 
mains inchoate. 

The "maximality" of a galaxy is often assessed via the 
st ellar-disk mass fraction, T^ lsk = V^ lsk /V c (equation 11 
in iPaper it ) . This quantity is traditionally measured at 
2.2h Rl the radius at which the circular speed peaks for 
razor-thin, radially ex ponential disks (see the generaliza- 
tion to oblate disks in IPaper Ij ) . In the idealized case of 
a two-component galaxy with a spherical DM halo and 
an exponential stellar disk, J^ lsk (i? = 2.2h R ) uniquely 
quantifies the influence of the DM on the mass budget 
at all radii for a given disk scale length, oblateness, and 
DM-halo density parameterization. This idealized case 
provides a useful fiducial model with which to compare 
observations, as we discuss below. Real galaxies deviate 
from the idealized case due to (1) the inclusion of other 
baryonic components with generally different mass distri- 
butions, such as bulges and gaseous disks, and (2) per- 
turbations of the stellar-disk mass-surface-density pro- 
file away from the nominal exponential, as inferred from 
surface-brightness variations. Thus, as explicitly asso- 
ciated with the stellar disk, J^ lsk (i? = 2.2h R ) has an 
intrinsic distribution for maximal disks: Sackett (1997) 
adopted sk (R = 2.2h R ) = 0.85 ± 0.10 as an appro- 
priate definition for maximal disks in galaxies of similar 
Hubble type to the Milky Way (Sb to Sc), also represen- 
tative of the DMS Phase-B sample. We directly compare 
this definition to our measurements in UGC 463; how- 
ever, we note that the literature studies upon which this 
definition was based do not remove the molecular gas 
component from the total disk mass distribution as we 
do for UGC 463. 

As roug hly synonymous throughout the discussion 
bv lvAS86L the "maximum-disk" hypothesis could also 
be termed the "maximum-baryon" or "minimum-dark- 
matter" hypothesis. In this respect, it is also useful to 
calculate the baryonic mass fraction, T h = As 
implemented byjyAS86, J 7 ^ « J\^ lsk at least in the sense 
that their mass-to-light ratios incorporated all mass dis- 
tributed similarly to the luminous disk and they assumed 
the stellar disk was by far the most massive baryonic 
component. Our detailed accounting of multiple bary- 

26 This should not be confuse d with the total baryon fraction, 
^"bar = M^JM^ n , discussed in IPaper O . 



onic components in UGC 463 with different mass distri- 
butions, particularly with regard to the massive gas disk, 
means that (1) F h ± J"f sk and (2) quoting LFf sk {2.2h R ) 
has a more limited bearing on the relative influence of 
the baryonic and DM mass on the total mass budget 
than described above. Therefore, it is useful to consider 
the radial functions J"^ lsk (i?) and J-b(-R) and to define 
multiple fiducial radii for J-"b based on the expectation 
of that each baryonic component can dominate the mass 
budget in distinct radial regimes. 

One expects J-b to decrease with radius as DM in- 
creasingly dominates the mass budget. The "maximum- 
disk" hypothesis effectively states J-^, ~ 1 at small ra- 
dius, regardless of whether or not the disk or bulge dom- 
inates the baryonic mass (see the imp lementation of th e 
maximum- disk" hypothesis by, e.g.. lKeml[T986L [1981 . 
The application of J-~b ~ 1 at small radius to Sb-Sc galax- 
ies has shown J r ^ lsk (2.2/i/j) = 0.85, which may be con- 
sidered a lower limit for j-y, at these radii. Finally, it 
is expected that Fb — > at large radius because rota- 
tion curves have been shown to remain nearly constant 
up to the radial extent of H I disks. The recent study 
by iDutton et al.l (|2011bf ) has shown that this expecta- 
tion for J-b(-R) is likely a limited picture of the range in 
galaxy properties. By com bining strong lensing anal- 
ysis and stellar kinematics, IDutton et al.l have shown 
T h (R e ) = 0.99 and T h {2.2h R ) = 0.67 for the late-type 
lens galaxy SDSS J2141-0001, which is at a redshift of 
0.14; R e is the effective (half-light) radius of the bulge 
(Section l2.2.2p . Thus, while adhering to the fundamental 
tenant of the "maximum-disk" hypothesis, SDSS J2141- 
0001 exhibits a disk that is relatively less massive than 
the disks of local Sb-Sc galaxies, such that it would be 
co nsidered "subm aximal" under the definition proposed 
bv lSackettl (|1997h . 

For UGC 463, we assess the ratio of the dark and bary- 
onic mass components via the mass fractions J-"^ lsk and 
J-h in Section IB. 3. 11 We quote Tb at R e — 0.5 kpc and 
2.2h R = 7.8 kpc as to compare with both the expecta- 
tions fro m the "maximum-disk" hy pothesis and the re- 
sults from IDutton et aTl ^Ollbf l Fl These radii are close 
to the radii at which the circular speeds of the stellar 
components are maximized (Figure I14|) . We also quote 

27 In galaxies with extended gas disks, one might also quote 
J-j, at a radius where t he g as disk dominates the baryonic circular 
speed; however, Figure 1141 shows that there is no such radius for 
UGC 463 within the radial regime of our measurements. 
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jrdisk^2_2ft, i? ) for a direct comp arison wi t h the definition 
of a maximal disk proposed by Sa ckettl (|1997l) . Further 
assessments are made using the ratio of the mid-plane 
volume densities (pdm/ Pb)z=o in Sec tion 16.3.21 and the 
enclosed-mass budget in Section l6.3.3l In all of these sec- 
tions, we compare our measurements for UGC 463 with 
an idealized, maximal-disk galaxy; this model provides 
the equivalent quantities for a maximal disk that has 
been embedded in the NFW and pseudo-isothermal ha- 
los fitted to our measurements of Vdm (Section l6.2p . This 
fiducial maximal disk has the scale length and oblateness 
as measured for UGC 463, is purely exponential in both 
R and z and has Jdisk = 0-85 at R — 2.2hn. Here, our 
notation Jdisk, as opposed to J-"* Isk , purposely leaves the 
nature of the matter in the disk undefined, as done by 
vAS86 in their implementation of the "maximum-disk" 
hypothesis. The combination of this fiducial maximal 
disk and our fitted DM halos produces a substantially 
higher circular speed than we measure for UGC 463, 
and the shape of the rot ation cu rve does not adhere to 
the disk-halo conspiracy (lvAS86D : we have shown in Scc- 
tion l6.1l that application of the "maximum-disk" hypoth- 
esis to UGC 463 essentially results in a marginal DM-halo 

mass within the radial regime pro bed by our data . 

Before continuing, we note that iBottemal (|1993l 119971) 
used kinematic measurements similar to our own and 
found 7 r disk(2.2/i^) = 0.63 ± 0.10 for a sample of 12 late- 
type galaxies; Bottema did not differentiate between stel- 
lar and gas components in this decomposition, hence our 
use of .Fdisk- In detail, Bottema's assumptions are not ex- 
actly the same as our own; however, we can cast UGC 463 
in terms of his measurement to find J r di S k(2.2/i j R) = 0.6, 
which is compatible with these previous results. We also 
not e that, using planetary nebu lae as kinematic trac- 
ers, iHerrmann fc Ciardullol (|2009t ) reported submaximal 
disks for four of the five galaxies they studied. 

6.3.1. Stellar-Disk and Baryonic Mass Fractions, J\: 1S 

and J~b 

Figure [16] provides J-~* lsk and J\, as a function of radius 
assuming a constant T* t K(R). Individual data points use 
our direct measurements of V c , whereas the gray lines 
calculate U c 2 = V h 2 + V£ M using Vb from Figure [T4] and 
Vdm from the NFW and pseudo-isothermal halos fitted in 
Figure [T5l Both T^ lsk (R) and Jb(i?) are effectively con- 
stant between 1.0 < R/fiR < 3.5; the roughly constant 
J\> ~ J 7 * lsk at these radii reflects the similarity between 
the radial distribution of the gas and stars. Figure [16] 
also provides the expected .Fdisk for the fiducial maxi- 
mum disk, which is the same for the NFW and pseudo- 
isothermal halos at R > 0.5hn. The shape of J^ lsk (i?) is 
very similar to that calculated for a maximal disk, which 
is essentially a statement that [i' K is very close to an expo- 
nential after subtracting the central light concentration; 
however, the normalization is very different. At small 
radius, T\> and J 7 * lsk diverge largely due to the exclusion 
of the central mass concentration in the calculation of 
the latter. 

TablelEgives T h {R e ), T h {2.2h R ), and sk (2.2h R ) for 
both the constant and variable T* ,k{R) assumptions. 
Our measurement of J r b(i? e ) has significant uncertainty 
due to the error in measurement of V c at this radius 
(largely due to centering errors), the uncertainty in the 
baryonic rotation speed, and the difference in the ex- 
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Figure 16. Mass fractions (T = V/V c ) of the stellar disk (J"^ isk ; 
filled points) and all baryonic matter (F^; open circles) as a func- 
tion of radius assuming a constant T» if(il). The discrete mea- 
surements use the measured values of V c ; random errors are shown 
in black for both quantities, whereas systematic errors are shown 
in gray only for measurements of T\, . Gray lines show the result of 
calculating V c 2 = V^ M + V? using the NFW (solid gray lines) and 
pse udo-isothermal (dashed gray lines) DM halo models from Figure 
1151 Black lines show J^isk for a fiducial maximal disk (J-disk = 0.85 
at R = 2.2hn) that has been embedded in the fitted DM halos (see 
text). 

trapolation when based on either the NFW or pseudo- 
isothermal halo. To the contrary, both and .F* lsk at 
2.2hn are relatively well constrain ed. Following the def- 
inition proposed by ISackettl ()1997t ) and assuming a con- 
stant T*_k(R), these quantities demonstrate that UGC 
463 has a substantially submaximal disk by a factor of 
~ (0.85/0.46) 2 w 3.4 in mass; this reduces to 1.9 if one 
instead defines J 7 b(2.2hji) = 0.85 as a maximal disk. 
These factors are consistent with our previous discus- 
sion in Sections 15.6.21 and 16.11 regarding the comparison 
of ou r dyn amical measurements with the SPS predictions 
from lB03l Adopting a variable T* .k{R) results in a mass 
profile of the disk that yields a peak rotation speed at or 
beyond the limit of our calculation (i? = 15 kpc), whereas 
the constant T^^iR) disk has a peak rotation at 2.7h R . 
However, in both cases the stellar (and baryonic) disk 
remains submaximal due to the relatively constant value 
£ jrdisk ^ an( ^ jr b ) a ^ ji > j lR Although it is possible that 
the baryonic mass is close to satisfying the Fh{R e ) ~ 1 
given the uncertainty, our data suggest that the circular 
speed may have substantial contributions from DM even 
within the bulge region. 

6.3.2. Mid-Plane Volume- Density Ratio, (pdm/ pt>)z=o 

Our mass-surface-density and pdm measurements from 
Figures Q3] and [T51 respectively, provide the mid-plane 
volume-density ratio between the dark and baryonic mat- 
ter, (pdm/ Pb)z=o- We assume that the stratification of 
the gaseous and stellar disks combine to produce an 
exponential vertical density distribution with a single 
scale height h z and that the bulge is spherical. Thus, 
the baryonic volume density at the disk mid-plane is 
(Pb) 2 =o = (S st -r-E gas )/2/i z - r -/9buige- If the gas is more con- 
fined to the plane, this calculation results in upper limits. 
Figure [17] provides the mid-plane volume-density ratio 
when assuming a constant T +i if(i?), where the discrete 
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Figure 17. Mid-plane density ratio (ptita / Ph) z—o °f au baryonic 
matter as a function of radius assuming a constant k(R). The 
calculation of (pb)z-o is discussed in the text. Data points use 
the Pdm measurements from Figure 1151 whereas gray lines use the 
NFW (gray solid lines) and pseudo-isothermal (gray dashed lines) 
parameterizations of Pdm- Only random errors are shown; upper 
limit arrows are used when the error is larger than the measured 
value. The results for the fiducial maximal disk are plotted as black 
lines for both the NFW (solid black lines) and pseudo-isothermal 
(dashed black lines) DM halo parameterizations. 



measurements use the non-parametric calculation of pum 
and the gray lines use the NFW and pseudo-isothermal 
sphere parameterizations. We also plot the expectation 
for t he fiducial maximal disk discussed above. For refer- 
ence, [Bienayml[eFaI] (|2006() find (pr>M/pb)z=o = 0.14 in 
the solar neighborhood assuming a spherical halo. 

At R 3> z, the change in pdm with z is much smaller 
than the change in p^, such that (pdm/Aj) oc exp(\z\/h z ). 
Averaging the data at 1 < R/hji < 3 in Figure H71 we 
find that (pdm/ Pb)z=o ~ 0.2, which is a factor of 5 larger 
than the expectation for a maximal disk (only 40% larger 
than the Milky- Way value). The mass volume density 
of UGC 463 is, therefore, dominated by dark matter at 
\z\ > l.6h z at R> Iir. This result is particularly impor- 
tant for our understanding of out-of-plane motions in the 
disk of UGC 463: The derivation of S dy „ = a 2 z /TrkGh z 
assumes an isolated, plane-parallel, infinite disk. Devi- 
ations from these assumptions, such as embedding the 
disk in a very massive DM halo, introduces s ystematic 
errors in the calculation, as briefly discussed in lPaper III . 
The result for (pdm/ Pb)z=o in UGC 463 suggests that 
s uch effects may b e significant for this galaxy. 

Bottcma ([1993D has discussed the influence of a mas- 
si ve DM halo on a z in disk stars, continuing the work 
of IBahcaU (|1984f) . These authors find that a z should be 
inflated relative to an isolated disk when embedded in 
a massive, spherical halo; the degree of the inflation is 
proportional to (pbm/ Pb)z=0j as shown in Figure 15 from 
iBottemal (I1993J). For th e fiduc ial maximal disk shown in 
our Figure [T71 Bot temal (|1993f) would predict a less than 
5% increase in the a z over an isolated disk, whereas our 
measurements for UGC 463 from Figure [T7] suggest a z 
could be increased by 10 — 15%. This means that our 
calculation of Sdyn °c <7» could overestimate the mass of 
the disk by 20 - 30%. 

Ideally, one would calculate Edyn by first assuming an 



isolated disk and then converging to a solution that in- 
corporated the effects of the DM halo. However, we have 
not done so here for UGC 463 because (1) the random er- 
ror in our isolated-disk measurements of £d y n are of the 
same order as this systematic correction; (2) there is sub- 
stantial error in our measurement of (pDM/pb)z=o ~ 0.2, 
even if we adopt the parameterized solutions for pdm! 
and (3) there are equally unknown competing system- 
atics that work in the opposite direction, such as the 
inclusion of a massive, razor-thin gas disk. The contin- 
ued study of the effects of a massive halo on the ve- 
locity dispersion, as oppo s ed to just the rotation curve 
(|Havashi fc Navarro! [20061 iWidrowi [2008) . is worthy of 
a dedicated effort. However, given ambiguities regarding 
the three-dimensional structure of galaxies and the verti- 
cal stratification of disks (see discussion in Section [5.6.ip . 
a detailed understanding of the influence of the DM halo 
on the disk is complicated. Here, we simply note that 
the impact of a relatively massive DM halo in UGC 463 
works to further lower the maximality of an already sub- 
maximal disk. 

6.3.3. Enclosed-Mass Budget 

As discussed in IPaper it . the total baryon fraction 
-^bar = M h ot / M^y^ is ill-defined; however, the surface 
density and Vdm measurements from Figure [141 allow for 
a robust calculation of the enclosed-mass budget to a fi- 
nite radius, assuming the halo is spherical. The resulting 
mass budget is presented in Table [7] and the mass growth 
curves are shown in Figure [T5J We end the calculation 
at R = 15 kpc (R = i.2h R ; ~ i? 2 oo/10 for our fitted 
NFW halo), well within the limiting radius of our pJ K 
measurements (Figure but extrapolating beyond our 
dynamical data. We calculate the DM mass at 15 kpc us- 
ing the fitted halo parameterizations and the percentage 
errors from the measured data. The mass growth curve of 
the fiducial maximal disk discussed above is overplotted 
in Figure [18] for reference. Although Table [7] only pro- 
vides the results when assuming a constant T*^ (R), the 
results assuming a variable T* k(R) are insignificantly 
different. 

Figure [IS] shows that the integrated baryonic mass of 
UGC 463 is, at most, equal to the integrated DM mass 
at R ~ Iir. Beyond this radius, the integrated DM mass 
quickly begins to dominate, such that Mb/MnM ~ 0.3 
at 15 kpc. To the contrary, the integrated mass of the 
fiducial maximal disk always dominates over the DM 
mass; the minimum mass ratio TWdisk/A^DM ~ 1-4 is 
at R = 15 kpc. 

7. SUMMARY 

In this paper, we have presented a detailed case study 
of the dynamics and implied mass budget of the low- 
inclination, SABc galaxy UGC 463. We find the galaxy 
to be dominated by DM at nearly all radii, a classifica- 
tion usually reserved for low-surface-brightnes s galaxies 
where as UGC 463 is ~ 1 magnitude abo ve the | Freemanl 
(1970) mean central surface brightness ([Paper J ). The 
stellar disk of UGC 463 is submaximal by a factor of 
> 3 in ma ss. Submaximal disks have been b oth indi- 
rectly Ce.g. iCourteau fe Rixljl999f) or directly ([Bottemal 
ll993tlHerrmann fc Ciardullol l2009ft measured in the past, 
and our results are consistent with these measurements. 



Westfall et al. 



30 



R/h R 
0.5 1 




R (kpc) 

Figure 18. Mass growth curves for baryonic mass (solid black 
line), NFW (gray solid line) and pseudo-isothermal (gray dashed 
line) DM halo parametcrizations, and the fiducial maximal disk 
(dotted line). The baryonic mass profile assumes a constant 
Y*,k(R)', the dark- and light-gray regions represent the random 
and systematic e rror s, respectively. Individual measurements of 
Vdm from Figure IT5l are converted to masses assuming a spherical 
halo and plotted as black triangles; black and gray bars represent 
the random and systematic errors, respectively. 

At present, this general result should be unaffected by 
systematic error, despite the albeit large number of as- 
sumptions. A summary of our analysis follows. 

Our analysis of UGC 463 draws from nearly all of our 
survey data for this galaxy as collectively described in 
Section[2 We derive a distance of D = 59.67±0.01±4.15 
Mpc using the flow-corrected systemic velocity and Hq = 
73±5 kms -1 Mpc -1 . We derive optical and NIR surface 
brightness profiles from, respectively, archival SDSS and 
2MASS data; we use the latter to calculate a total K- 
band magnitude of rax = 9-32 ± 0.02. Our photometry 
is corrected for Galactic extinction; internal extinction 
and dust emission are marginal considerations for our 
results. We also correct for instrumental smoothing of 
the surface brightness profile at small radius. Ionized-gas 
and stellar kinematics are derived using SparsePak and 
PPak IFS. Ionized-gas kinematics combine results from 
atomic emission lines near Ha and the [O iii]A5007 line; 
kinematics are measured using single or double Gaussian 
line fits (as in I Andersen et aLll2008[ ). Stellar kinematics 
are derived from absorption lines with rest wave lengths 
between ~ 492 — 522nm using DCS ([Paper IIJ ) and a 
single Kl III template star (HD 167042 for SparsePak 
and HD 162555 for PPak): template mismatch is < 5% 
(|Paper III ). All kinematics are corrected for instrumen- 
tal broadening. We measure H I mass surface densities 
and the H I velocity field using 21cm aperture synthesis 
imaging from the VLA. Finally, we estimate the H2 mass 
surface density by combining 24/im Spitzer imaging and a 
hipm /IgnA V ca libration based on our reanalysis of data 
tabulated bv lL08l Errors in our measurements of Sh 2 in- 
corporate the systematic error in this calibration and the 
error in our adopted value of Xco = (2-7 ± 0.9) x 10 20 
cm -2 (K kms -1 ) -1 . 

In Section 13 we determine the detailed on-sky geomet- 
ric projection of the disk of UGC 463, including extensive 
tests of the inclination. These efforts are important to 



subsequent analysis of the measured kinematics due to 
the substantial influence of inclinatio n errors o n the fun- 
damental calculations in this study ([Paper IIP . Inclina- 
tions are measured both kinematically — - using an algo- 
rithm explained in lAndersen et ail ([20081 ) that assumes 
circular motion and a single, copla nar disk — and by 
inverting the TF relations from IVOTI . Kinematic inclina- 
tions are most consistent between all ionized-gas and stel- 
lar tracers when adopting a d Vlos / di weighting scheme, 
as explained in Appendix [Bj We find «kin = 25? 1 ± 2?5 
and «tf = 29° ±2°; a combined measurement of i — 
27° ± 2° represents our best estimate for the inclination 
and is used throughout all subsequent analysis. We de- 
rive the dynamical center of each SparsePak observation 
and, for the PPak data, we affix the dynamical center 
to the morphological center determined from a recon- 
structed continuum map; our two-dimensional maps in 
Figure [T] demonstrate that the dynamical and morpho- 
logical centers of UGC 463 are identical to within the 
errors of our measurement (~ 1"). 

Using the derived geometry, we create azimuthally 
averaged kinematics in Section [4] To combine kine- 
matic measurements using different instruments and dif- 
ferent wavelength regimes, we apply beam-smearing cor- 
rections to both our ionized-gas and stellar kinematics; 
be am-smearing correc tions to the H 1 data are described 
by iMartinsson (| 2 1 If ) . The beam-smearing corrections 
employ model surface-brightness, velocity, and velocity- 
dispersion distributions to create a synthetic dataset that 
is compared to our observations; the corrections are small 
(less than a few percent) except for data near the dy- 
namical center. Figure [7] demonstrates that both the 
ionized-gas and stellar kinematics measured separately 
by SparsePak and PPak are very well matched. We pro- 
vide a cursory assessment of the kinematic axisymme- 
try of the rotation curves and velocity dispersion pro- 
files by overlaying 180° azimuthally averaged kinemat- 
ics for the receding and approaching sides in Figure [7J 
with respect to the errors, only moderate differences are 
present. Therefore, we impose axisymmetry by measur- 
ing the gas (ionized+neutral) and stellar velocity and 
velocity dispersions only as a function of in-plane galaxy 
radius. We correct the gas rotation curve to the circu- 
lar speed using me a surem ents of a gas and E gas following 
iDalcanton k, Stilpl ([2010D to produce the axisymmetric 
radial profiles in Figure [5] 

Based on our collection of azimuthally averaged prop- 
erties, we determine physical properties of the disk of 
UGC 463 in Section [5] We find that the stellar ve- 
locity dispersion profile, ~aZ(R), has only minor devia- 
tions from a pure exponential; the fitted e-folding length 
(her ~ 2.6hn) suggests either Td yn or h z increases by 
a factor of ~ 2.3 over the radial range of our data, if 
the other quantity is radially invariant. Using the mea- 
sured circular speed, stellar rotation curve, and stellar 
velocity dispersion profile, we calculate the observable 
function Sad(R)- This function is directly related to the 
shape of the SVE such that we find a = 0.48 ± 0.09 
and f3 = 1.04 ± 0.22, when assumed to be constant over 
the entire disk. Therefore, we find the conversion factors 
<Jz/(<?*)e = 0.76±0.09 and a R /{a*) e = 1.59±0.17, which 
we use in our calculations of Edyn ano - ^ ne disk stability, 
Q. Our calcul ations of Edyn (Figure [10]) use equation 9 
from IPaper 111 , which assumes the oblateness from their 
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equation 1. We find Edyn is well fit by 3£ gas , from which 
it follows that £» ~ 2£ gas . Using £ gas , c gas , £*, and or, 
we calculate individual and multi-com ponent s t ability 
coef ficients followin g from, respectively, iToomrd ()1964[ ) 



and lRafikovl (|2001D . We find the disk to be globally sta- 
ble, with the multi-component stability asymptotically 
decreasing with radius to a value of Qr ~ 2. Based on 
combini ng stability arguments with swing-amp lification 
theory (jToomrel fl98lt lAthanassoula et aLlll987t) , disks 
of a fixed rotation curve should exhibit higher spiral-arm 
multiplicity when the disk mass is decreased. This expec- 
tation is qualitatively consistent with our measurements 
of a submaximal disk and t hree- arm mu l tiplici ty in UGC 
463; however, Paper V and Martinsson (2011) find sub- 
maximal disks regardless of spiral-arm multiplicity. 
We calculate T£g using £* and fi' K to find (T^) = 

0.22 ± 0.091^ M Q /£% at R > 2 kpc. Our measure- 
ments of (T* 1 *) are systematically lower than SPS model 
predictions (T^) from IB03I by a factor of ~ 3.6 and 
larger than the SPS modeling of Z09 by a factor of ~ 1.8. 
Based on an MC sampling of probability distributions 
assigned to each quantity in the calculation, we gener- 
ate a composite (random+systematic error) probability 
distribution for T^f (R) and (T^f >. We find that the 
Z09 prediction are within our 68% confidence interval for 
(T^ 1 *); in contrast, measurements of (T^ 1 ^) consistent 

with the[B03 prediction occur for less than 1 in 10 5 MC 



samples. IZ09I attribute the disparity between their Tf p ^ 
predictions and those from lB03l to different treatments of 
the star-formation history and TP-AGB phases of stellar 
evolution. We also find a factor of ~ 2 increase in T^ 1 ^ 
with radius, which is not predicted by the SPS modeling. 
This feature may reflect a true increase in the T* 1 * , a 
flaring of the stellar disk, or a change in the relative dy- 
namical influence of the halo, thick stellar disk, and/or 
razor-thin gas disk. 

We discuss the mass budget of UGC 463 out to 15 kpc 
(4.2 hp;) in Section [6] using a traditional rotation-curve 
mass decomposition, which benefits from our unique and 
direct measurement of T^ 1 ^. Our primary discussion as- 
sumes a constant T»j^(i?) = (T^ 1 *); however, we also 
briefly discuss results obtained by assuming a variable 
^*,k{R) = "^'^(i?). Our mass decomposition also as- 
sumes that h z is constant with radius and that the galaxy 
is composed of four separable potentials (halo, stellar 
bulge, stellar disk, and gas disk). The total mass budget 
obtained by assuming a constant T*,j<- (R) is provided in 
Table [3 assuming a variable T^^i?) produces results 
that are statistically identical. 

We calculate the circular speed of each baryonic com- 
ponent in Section 16.11 We find that a maximal disk 
may be produced by adopting the prediction from 

IB03I . which amounts to increasing the stellar (baryonic) 
mass by a factor of 3.6 (2.7) above our measurements of 
(Yf'sij). jj oweveri this result is effectively excluded by 



- *,K/ 

our simultaneous measurements of V c and o*. Using our 
dynamical measurements, we produce V^ M — U c 2 — U b 2 
and use these measurements to calculate pdm; both Vr>u 
and pbm assume a spherical halo. In Section 16.21 we 
fit Vdm with an NFW and pseudo-isothermal DM halo, 



and we find both to be statistically suitable descriptions 
of our measurements. The concentration of the NFW 
halo is consistent with expectations from DM-only sim- 
ulations, implying that the halo structure has been rela- 
tively unaffected by the collapsed baryons. Our measure- 
ments of pdm show a slope that may be steeper than both 
Pdm oc 1/R and the slope predicted by either DM- halo 
parameterization; however, this result is highly depen- 
dent on the error-prone assessment of the baryonic mass 
within the central kpc. Results for additional galaxies in 
our Phase-B sample are required to place better statis- 
tical constraints on the shape of DM halos of local disk 
galaxies. 

We discuss the dominant gravitational influence of DM 
over baryonic matter in UGC 463 in Section 16.31 as 
summarized in Tables [7] and |5] Assuming a constant 
Tf*,if(i?), we find the baryonic disk to be substantially 
submaximal with F h {2.2h R ) = 0.61±£'o9 toil- Consid- 
ering only the stellar disk and adopting J r f lsk (2 .2hp) = 
0.85 ± 0.10 as the definition of a maximal disk (|Sacketti 
1997), UGC 463 is submaximal by a factor of ~ 
(0.85/0.46) 2 = 3.4 in mass, consistent with our expec- 
tation based on the difference between (TjJ 1 ^) and the 
T^jf prediction from IB031 (Section 16.11) . We also com- 
pare our measurements of the baryonic component to a 
fiducial maximal disk (having .Fdisk(2.2/iR) = 0.85) re- 
sulting from embedding a purely exponential disk (with 
hp and h z as measured for UGC 463) in our fitted NFW 
and pseudo-isothermal DM halos. In the disk mid-plane, 
we find the ratio (/?DM/pb)z=o is a factor of five larger 
than expected by our fiducial maximal disk, which may 
lead to an overestimate of Sdyn- Additionally, we find 
that the enclosed-mass of the galaxy is dominated by 
DM at R > hu, whereas the fiducial maximal disk dom- 
inates the enclosed-mass budget at all radii sampled by 
our observations. Finally, assuming a constant Y*^k(R), 
we find a baryonic-to-DM mass ratio of 0.31 ± 0.81^ 
for the mass enclosed within the central 15 kpc. 

Although our results are for a single galaxy, an anal- 
ysis of 30 galaxies in our sample demonstrate that all 
of these disks ar e comparably submaximal (Paper V; 
iMartinssonl 120 lift . One can increase the maximality of 
these disks by changing the assumptions concerning the 
vertical mass distribution (quantified by the constant k) 
or the applied oblateness (q), as discussed in Sections I5~3l 
and l5.6.T1 however, k and q would have to take on values 
that are effectively excluded by empirical constraints on 
the structural pa rameters of d isk galaxies based on edge- 
on systems (e.g.. IKregell 120031) . The recalibration of T* 
based on these results has significant consequences for, 
e.g., our understanding of the baryonic mass of galaxies 
as a function of redshift and for the gravitational inter- 
play between baryonic and dark matter in the process of 
galaxy formation. 
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APPENDIX 

A. SURFACE-BRIGHTNESS CALIBRATION AND INTERPOLATION OF IFU DATA 

The spectral-continuum surface-brightness maps in the Ha and Mgi regions shown in Figure Q] are calibrated against 
SDSS imaging data, assuming that the dynamical center (Table 15]) is the same as the morphological center. 
The "model" flux in fiber / measured from the CCD data is 

C mod ,/ = dlog [0A(Z nb - Zccb) + log/ccD,/] A &h - S fi b, (Al) 

where Z represents a magnitude zero-point and Iqgdj = (Cccd,/ — SccdV^ccd represents the surface brightness in 
units of DN arcsec -2 — determined by the total flux (C) within an aperture (of area A) with a sky background (S). 
The fiber aperture area, Aqcb ~ ^-fib = ^D^/A, is known and the quantities related to the CCD image are measured 
or provided by the SDSS calibration. The fiber-continuum zero-point, Z^, and sky-level, Sab, are free para meters; 
iSfib is fiber independent, adjusting the nominal correction based on the average sky spectrum. Equation lAll is fit to 
our IFS data by minimizing (yJC^hj — ^C nw dj) 2 , where Cfib,/ is the mean flux across the full spectral range; data 
with erroneous negative flux are ignored. We also limit the radial region considered to avoid inflated errors were the 
sky subtraction of the IFS is particularly problematic due to variations in the sky flux as measured by the dedicated 
sky fibers. 

We use the SDSS g-band (1'.'5 seeing) and r-band (172 seeing) data to calibrate the Mg i-region and Ha-region 
IFS, respectively. For the PPak data, analy sis of guide-camera images taken throughout each exposure show that 
the average seeing was l'/7 ( M artinssonl 1201 II ) : therefore, we match this seeing by applying a Gaussian kernel with a 
FWHM of 078 to the <?-band image when fitting to the PPak data. No such seeing measurements are available for the 
SparsePak data. Although we have allowed the seeing to be a fitted parameter for SparsePak, seeing measurements 
were non-convergent in the sense that there appears to be no substantial difference with the inherent seeing of the 
SDSS images. This is not surprising given the image quality quartiles at WIYN and the large SparsePak fibers. 

The final calibration results are provided in Figure [19] for all data. We mark each panel by the pointing number and 
include the residual RMS within the fitting region. The RMS values are typically 0.2 magnitudes, and the PPak data 
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Figure 19. Photometric calibration of the IFS fiber continuum flux. Three panel groups are shown: the results for the four Ha 02Jan02 
SparsePak pointings, the results for the three He* 20Oct02 pointings, and the results for the two Mgl pointings — one from SparsePak and 
the other from PPak. Each group has two columns, overlaying the calibrated fiber fluxes (black) on the measured aperture photometry 
from the SDSS images (gray) to the left and plotting the residual, A/i = jUfib — MCCDi to the right. Radial regions not considered during 
the calibration process are shaded gray. The pointing number or instrument is provide in the upper-right corner of the left-column panels, 
and the RMS of the residuals are displayed in the upper-left corner of the right-column panels. 

have the smallest residual at 0.1 magnitudes. Fits using no additional seeing for the SparsePak data demonstrate good 
agreement with the direct-imaging data at small radii, implying that the systematic errors due to an inappropriate 
seeing match between the SDSS images and the SparsePak IFS are inconsequential. The results shown in Figure [TH] 
are used to produce calibrated fluxes that are interpolated and converted to surface brightness for Figure [T] 

The interpolation of the continuum fluxes and kinematics from our IFS is performed to fill interstitial fiber regions 
according to the following algorithm: Each kinematic measurement contributes to every pixel in the image, weighted 
by a two-dimensional Gaussian function centered on the fiber aperture with a FWHM equal to a multiple of the 
effective fiber diameter, -D n b, and inversely weighted by the measurement error; therefore, each image pixel is weighted 
both by the quality and proximity of the surrounding fiber measurements. Mathematically, the weight of each fiber / 
at pixel coordinate (J,k) is, therefore, 



w f (j,k) 



1 



■ exp 



(nAi b /2) 2 /ln2 



(A2) 



where e/ is the measurement error, (xj,yk) are the on-sky coordinates of pixel {j,k), the fiber center has on-sky 
coordinates (xf,yf), and n£>fib is the FWHM of the circular Gaussian in multiples of the fiber diameter. Calibrated 
continuum fluxes include no additional error weighting (e/ is constant for all /). The interpolated value at each 
pixel is then the weighted average, over all fibers, of the continuum value or kinematic measurement in question. 
For presentation purposes, we limit the interpolation to only those regions with a "significant" contribution to the 
interpolated map. For Figure [U we adopt the following representation of "significant:" 



J2 w f(i> k ) ^ 
/=i 



N X Ny . 



Nx-l Ny-1 N f 

12 12 12 w fti> k ^ 

k=0 f=l 



(A3) 
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Table 9 

Interpolation 
Parameters for Figure [T] 



Quantity 


n 


Wmm 




1.4 


0.5 


v Ka 


1.4 


0.1 


WVIgi 


1.2 


0.5 


V, 


1.8 


0.1 




1.8 


0.1 



where N x and N y are the pixel dimensions of the image. That is, the summed weight of all fibers to a given pixel 
must be greater than w m i n times the mean of all weights across the entire image. This scheme is not ideal given the 
dependence between the inclusion of an interpolated value and the arbitrary size of the image; however, it can provide 
reasonable results as demonstrated by Figure [T] Values of n and w m i n for each of the 5 interpolated maps in Figure Q] 
are provided in Table [9j 

B. OPTIMAL WEIGHTING SCHEME FOR KINEMATIC INCLINATION MEASUREMENTS 

In velocity-field modeling, parameter robustness and covariance can be, respectively, improved and mitigated by 
introducing data-weighting schemes. For example, iBegema 3 (|1989f ) introduced' a cosine weighting scheme, effectively 
weighting each datum by the derivative of the model LOS velocity with respect to the projected rotation velocity, 

w — dVhos /<9Kot° J = cos 6*, thereby limiting the covariance between Vrot and i. lAndersen k, Bersha dv (2003) mitigate 
this same covariance by instead fitting the projected rotation curve directly. Here, we consider an optimal weighting 
scheme for measuring the kinematic inclination of UGC 463 according to the approach described in Section f3. 1.1 1 

lAndersen fc Bershadyl (|2003l ) parameterized V^ t ° } (R) by a hyperbolic tangent function and adopted a velocity-error 
weighting scheme, producing a face-on TF relation that is well matched to samples of more inclined systems. Their 
error-weighting scheme combines, in quadrature, the measured velocity error with a "beam-smearing error" and a 
"stochastic error." The beam-smearing error is based on a fiber-by-fiber measurement of the variance in Ulos within 
the fiber aperture, thereby reducing the effect of patchy emission on the fit. The stochastic error is a single error 
assessed for every velocity measurement that reduces the influence of small-scale, incoherent non-circular motions on 
the fit. 

Here, we are primarily concerned with fitting inclination such that we test the success/failure of a given weighting 
function, as applied to UGC 463, via the correspondence/disparity of the inclinations determined from each of three 
tracers: (1) Ha from SparsePak, in a ; (2) [O III] from PPak, iro ml j an d (3) stars from PPak, i». We apply four 
weighting schemes: In addition to the error-based and cos 9 schemes described above, we include uniform weighting 
and a weighting scheme defined by w = 8Vlos /9i. The latter scheme affords those data with greater leverage on the 
fitted inclination a greater influence on the goodness-of-fit statistic. 

We omit data from consideration in our goodness-of-fit statistic in two steps. First, we omit all data with velocity 
errors that are greater than 15 kms -1 , eliminating 0.5%, 0%, and 1.0% of the SparsePak Ha, PPak [O in], and PPak 
stellar data, respectively. We note here that, after applying this omission, the mean velocity measurement errors 
are 1, 3, and 6 kms -1 for the SparsePak Ha, PPak [O in], and PPak stellar data, respectively. Second, we omit 
highly discrepant velocities by first fitting the data using the error weighting scheme and omitting data at high % 2 . 
Data are iteratively omitted while adjusting the model and the stochastic error until the error- weighted distribution 
of the data about the model follows a nominal Gaussian (Wcstfall 2009); in practice, no points are omitted with 
X 2 = {V — V m ) 2 / e(V) 2 < 10. This omission stategy does not bias our results toward, e.g., the initial guess parameters 
of the fit; instead, it serves to eliminate a ^-optimization bias driven by a few, highly discrepant measurements. 
Applying this procedure to our UGC 463 data eliminates an additional 11.6 %, 14. 0%, and 5.2% of the SparsePak 
Ha, PPak [O ill], and PPak stellar data, respectively. As described in Section 13.1.11 all velocity-field parameters are 
simultaneously fit to the remaining data with the dynamical center of the PPak data fixed to the morphological center. 
Identical data sets a re fit by each wei ghting scheme for each tracer, and we use a set of 500 bootstrap simulations 
(see Section 15.6.2 of iPress et al.|[2007h to determine the probability distribution for each fitted parameter. The error- 
weighted standard deviations of the velocities about the best-fitting models are typically 5, 6, and 7 kms -1 for the 
SparsePak Ha, PPak [O in], and PPak stellar data, respectively, with only small (5 — 10%) variations among results 
reached using the different weighting schemes. 

Figure [20] provides the bootstrap-based probability distributions and growth curves for the inclinations measured 
by each weighting scheme and each tracer; the best-fit inclination and the 68% confidence limits are tabulated in the 
Figure. For UGC 463, we find that weighting by dV^os/di produces inclination distributions that are the most similar 
between the three kinematic tracers. This is, therefore, the weighting scheme we have adopted in Section 13.1.11 to 
measure iki n - For all results except those based on the error- weighting scheme, the error- weighted mean of irq, i[oni]i 
and i* is compatible with our final adopted inclination of i — 27° ± 2° to better than the errors. 
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Figure 20. Inclination probability distributions determined from 500 bootstrap simulations for each tracer using four different weighting 
schemes. Left — Histogram of the returned inclination values. The fitted inclination for each tracer is tabulated in the upper-left corner of 
each panel. The line key for each histogram is shown in the lower-left panel. The SparsePak Ha histogram is shaded in light (transparent) 
gray and the PPak [O III] histogram is shaded in dark gray. Right — Growth curves of each histogram from the left panels. The line types 
are repeated from the left column. The dotted lines mark the median (0.5 growth) and the 68% confidence interval (0.16 — 0.84 growth). 



